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Results  of  an  experimental  study  on  the  dynamic  properties  of 
OFHC  copper  (1/2-hard)  are  presented.  Areas  studied  included 
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and  20,204  and  370®C) .  Longitudinal  and  shear  wave  velocities 
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profiles  showing  both  release  wave  behavior  and  attenuation 
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ABSTRACT 


Results  of  an  experimental  study  on  the  dynamic  properties  of 
OFHC  copper  (1/2-hard)  are  presented.  Areas  studied  included 
stress-strain-strain  rate  and  reverse  loading  behavior,  elastic 
constants,  equation  of  state,  compressive  and  release  wave  charac¬ 
teristics,  and  spall  fracture.  The  material  showed  approximately 
elastic-perfectly  plastic  behavior  under  uniaxial  stress  compres¬ 
sion  with  slight  rate  sensitivity  (0.005  to  950/sec  strain  rates 
and  20,204  and  370“C) .  Longitudinal  and  shear  wave  velocities 
at  20°C  were  4.757  and  2.247  mm/usec,  and  temperature  and  pres¬ 
sure  dependence  was  also  measured  and  various  elastic  constants 

were  calculated.  The  shock  wave  equation  of  state  (EOS)  was 

2 

measured  from  10  to  200  kbars  as  a„  =  1.5  +  344  u  +  143  u 

n  P  P 

The  EOS  was  also  calculated  from  the  elastic  constants  and  com¬ 
parisons  made  with  the  shock  wave  EOS  Compressive  wave  tests 
showed  a  ramped,  non-steady-state  elastic  precursor  with  very 
rapid  decay  (>  80%  in  2  mm) .  The  yield  level  from  compressive 
wave  profile  measurements  was  considerably  lower  than  that 
determined  from  uniaxial  stress  and  EOS  tests.  Shock  wave 
profiles  showing  both  release  wave  behavior  and  attenuation 
were  obtained.  The  impact  velocity  required  for  spall  fracture 
was  found  to  increase  with  decreasing  pulse  width. 
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FOREWORD  -  EXPERIMENTAL  DATA  VOLUMES 

The  prediction  of  reentry  vehicle  response  to  impulsive  loading 
resulting  from  energy  deposition  has  been  studied  extensively 
during  the  last  decade.  Analytical  models  and  computer  routines 
have  been  developed  to  assess  the  vulnerability  of  such  vehicles 
to  the  initial  loading  phase,  in  which  stress  waves  are  gener¬ 
ated  and  propagated  through  the  structure,  and  to  subsequent 
phases,  where  elastic  vibration,  plastic  deformation  and/or  frac¬ 
ture  may  occur. 

A  Defense  Atomic  Support  Agency  (DASA)  program  was  initiated 
with  the  objective  of  improving  computer  codes  for  the  predic¬ 
tion  of  damage  induced  by  X-rays  (PREDIX) .  The  PREDIX  metals 
program  has  combined  the  efforts  of  several  contractors  under 
the  direction  of  DASA.  The  primary  contributions  of  each  con¬ 
tractor  can  be  summarized  as  follows: 

General  Motors  Corporation — Measurement  of  material 
properties,  including  stress-strain-strain  rate  be¬ 
havior,  elastic  constants,  equations  of  state,  shock 
wave  profiles  and  spall  fracture. 

Effects  Technology  Inc. — Spall  fracture  tests  and 
analysis  (exploding  foil)  and  underground  testing. 

KMS  Technology  Center — Constitutive  modeling. 

Physics  International  Company — Energy  deposition  and 
spall  fracture  studies  (electron  beam) . 

Systems,  Science  and  Software — Energy  deposition  and 
front  surface  phenomena,  constitutive  modeling  and 
code  development,  and  underground  testing. 
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INTRODUCTION 


This  report  is  one  of  a  series  of  six  giving  results  of  a  study 
of  dyneunic  properties  of  materials.  The  principle  objectives 
of  this  study  were: 

•  Provide  measurements  of  dynamic  properties  of 
materials  to  be  used  as  inputs  to  model  develop¬ 
ment  for  shock  wave  propagation  and  dynamic 
fracture . 

•  Provide  data  to  establish  accuracy  of  computer 
code  predictions. 

•  Provide  physical  interpretation  of  experimentally 
observed  material  response  to  dynamic  test  condi¬ 
tions  . 

The  primary  areas  studied  included:  (1)  Stress-strain-strain 
rate  behavior  in  uniaxial  stress,  including  Bauschinger  effect; 
(2)  Elastic  constants,  including  temperature  and  pressure 
dependence;  (3)  Hugoniot  equation  of  state;  (4)  Comprcjssive 
and  release  wave  profiles,  including  elastic  precursor  decay 
and  v:ave  attenuation;  and  (5)  Spall  fracture  studies,  including 
pulse  width  dependence.  The  six  reports  in  this  series  are: 


Vol. 

I: 

Summary  of  Results,  DASA  2501-1  (AD712847) 

Vol. 

II: 

Experimental  Methods  and  Techniques, 

DASA  2501-2 

Vol. 

III: 

6061-T6  Aluminum,  DASA  2501-3 

Vol. 

IV: 

Axpha  Titanium,  DASA  2501-4 

Vol. 

V: 

OFHC  Copper,  DASA  2501-5 

Vol. 

VI: 

Tantalum,  DASA  2501-6 
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The  experimental  technif^aes  utilized  in  this  study  are  reviewed 
in  Vol.  II  (DASA  2501-2) .  Briefly,  the  principal  techniques 
applied  were: 

Stress-strain  Studies — A  laboratory-type  universal 
testing  machine  was  used  for  rates  <  0.1/sec,  a 

_3 

medium  strain  rate  machine  for  rates  of  10  /sec 
2 

to  10  /sec,  and  a  Hopkinson  bar  device  for  rates 
2  3 

of  about  10  to  5  X  10  /sec. 

Elastic  Constants — The  pulse  superposition  method  was 
used  for  measuring  longitudinal  and  shear  wave  velo¬ 
cities  as  functions  of  temperature  and  hydrostatic 
pressure . 

Equation  of  State — Gun-launched,  flat-plate  impact 
techniques  were  used  to  generate  uniaxial  strain 
conditions.  Velocities  up  to  0.6  mm/ysec  were 
achieved  with  63.5  mm  and  102  mm  single-stage  com¬ 
pressed-gas  guns,  and  velocities  of  0.3  to  2  mm/ysec 
were  achieved  with  a  64  mm  single-stage,  powder  gun. 
Hugoniot  data  were  obtained  with  x-cut  quartz  gage 
and  streak  camera  techniques. 

Wave  Propagation  and  Spall  Fracture — Flat-plate  impact 
techniques  were  also  used.  Wave  profiles  were  measured 
with  x-cut  quartz  gages  and  a  laser  velocity  interfero¬ 
meter. 


A  literature  search  on  dynamic  properties  of  materials  was 
carried  out  under  the  PREDIX  program  and  a  listing  of  recent 
publications  (including  abstracts)  on  dynamic  properties  of 
materials  is  given  in  Reference  1. 
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SECTION  I 

MATERIAL  PROPERTIES 


Copper  has  a  face-centered-cubic  structure  and  the  material 

on  which  tests  discussed  in  this  report  were  carried  out  is 

* 

an  oxygen-free,  high-conductivity  (OFHC  )  grade  certified  to 
be  99.99%  minimum  copper.  Chemical  analysis  identified  no 
other  elements  present  in  amounts  greater  than  0.0005%. 

This  copper  is  not  heat-treatable  and  different  strength  levels 
are  achieved  by  cold-working.  The  material  tested  was  purchased 
as  63.5  mm  diameter  bar  stock,  which  was  found  to  be  preferrable 
to  plate  stock  since  it  puts  essentially  no  limit  on  specimen 
thickness  and  gives  little  scrap  while  ensuring  that  all  speci¬ 
mens  are  from  the  same  parent  material. 

The  copper  bar  stock  was  nominally  1/2-hard,  which  corresponds 
to  a  handbook  tensile  yield  strength  of  2.5  kbar  at  a  "static" 
test  rate.  Measured  hardness  was  78Rp  (equivalent  to  74  Bri- 
nell  @  500  kg  load) .  A  photomicrograph  of  a  representative 
portion  of  the  bar  is  shown  in  Figure  1.  The  bar  axis  or  cen¬ 
terline  is  vertical  in  this  picture  and  the  width  shown  is 
4  mm  out  of  a  total  of  63.5  mm.  There  is  some  grain  elonga¬ 
tion  along  the  axis  with  about  a  factor  of  three  variation  in 
grain  size,  the  average  being  approximately  ASTM  5  (64  urn) . 

For  results  presented  in  this  report,  the  test  or  wave  propaga¬ 
tion  direction  was  axial  (i.e.,  vertical  in  Figure  1). 

The  measured  density  (p  )  of  the  copper  was  8.917  g/cc 
3  ° 

(0.322  Ib/in.  ).  For  use  in  equation  of  state  and  energy 


Trade  Mark,  American  Metal  Climax,  Inc. 
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Figure  1  OFHC  Copper,  63.5  mm  Diameter  Bar  Stock 
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deposition  calculations,  several  physical  constants  were  compiled 

(2-4)  * 

from  the  literature  and  are  listed  below: 

Volume  coefficient  of  thermal  expansion,  3: 

51.5  X  lo"®  0.006  X  lo"®  X  T  /“C 

Specific  heat,  Cp: 

0.0915  +  0.00003  X  T  cal/g 

Melting  temperature:  1080*0 

Boiling  temperature:  2595*C 

Latent  heat  of  fusion:  49  cal/g 

Latent  heat  of  vaporization:  1146  cal/g 
(T  in  “O 

A  portion  of  the  copper  bar  was  annealed  for  use  in  some  s tress - 
strain  and  compressive  wave  tests.  Grain  structure  of  the  an¬ 
nealed  material  is  compared  with  as-received  stock  in  Figure  2. 
The  annealing  treatment  was  1  hour  at  500*C  in  air,  followed  by 
air  cooling.  Hardness  after  anneali ng  and  machining  away  of 
the  oxidized  suri.ace  layer  was  35  Rp. 


A  very  thorough  compilation  of  thermodynamic  and  elastic 
constants  for  the  elements  is  given  by  K.  A.  Gschneidner 
in  Solid  State  Physics,  Vol.  16,  1964.  The  constants 
listed  above  generally  differ  slightly  from  those  given 
by  Gschneidner  due  to  consideration  of  more  recent  data. 
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Figure  2  As-Received  and  Annealed  OFHC  Copper 
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SECTION  II 

STRESS-STRAIN  STUDIES 


Shock  wave  and  structural  response  code  input  requirements  in¬ 
clude  constitutive  equations  relating  stress  to  strain  and  to 
strain  rate.  To  obtain  these/  measurements  are  made  of  uniaxial 
con^iressive  stress  vs.  strain  at  various  strain  rates  and  tem¬ 
peratures.  These  tests  give  yield  and  flow  stress  behavior, 
strain  rate  sensitivity  and  work  hardening  characteristics. 
Although  it  is  recognized  that  heating  rate  and  time-at- tem¬ 
perature,  as  well  as  final  temperature,  are  important  in  assess¬ 
ing  material  properties  under  short-time  (<  10  ^  sec)  energy- 
deposition  conditions,  the  stress-strain  data  reported  here  are 
for  relatively  long  heating  times  only.  Bauschinger  effect 
tests  were  conducted  to  provide  information  on  unloading  and 
SLibsequent  yield. 


UNIAXIAL  STRESS  TESTS 

Results  of  conpressive  stress-strain  tests  on  copper  are  shown 
in  Figure  3  for  strain  rates  of  0.005/sec  to  SSO/sec,  with  each 
curve  the  average  of  several  tests.  The  as-received  copper  was 
tested  at  three  temperatures  (20,  204  and  370*C) ,  while  a  few 
tests  were  performed  on  annealed  material  at  20*C.  The  heating 
rate  for  the  elevated  temperature  tests  was  'v  0.5  *C/sec,  and 
the  reduction  in  yield  and  flow  stress  observed  at  elevated 
teiqperature  may  be  a  function  of  heating  rate,  as  has  been  ob¬ 
served  for  6061-T6  aluminum. 
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Figure  3  Conpressive  Stress  vs  Strain,  Copper 
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The  strain  rates  given  are  average  plastic  strain  rates,  with 
strain  rate  variation  during  the  test  less  than  a  factor  of 
two  in  the  plastic  region.  The  effect  of  this  variation  has 
been  shown to  be  insignificant  when  considering  six  or 
seven  orders  of  magnitude  range  in  strain  rate.  Accuracy  of 
the  results  is  ±  3%  for  low  and  intermediate  rates  and  ±  10% 
for  the  high  rate.  Repeatability  was  better,  however,  and 
scatter  due  to  material  variability  and  technique  error  was 
less  than  2%  for  low  and  intermediate  rates,  and  less  than 
4%  for  the  high  rate  above  2%  strain. 


The  stress-strain  curves  in  Figure  3  for  as-received  (1/2-hard) 
copper  show  an  elastic-plastic  behavior  with  very  littZ  /ork 
hardening,  at  all  strain  rates  and  tenperatures  studied.  The 
emnealed  copper  shows  an  order  of  magnitude  decrease  in  yield 
strength  but  a  pronounced  increase  in  work  hardening.  Work 
hardening  exponents,  N,  were  calculated  from  the  compressive 
stress-strain  curves  by  applying  the  following  relation  to  the 
plastic  region: 


where  a  and  e  are  stress  and  strain  at  a  given  point 
in  che  plastic  i 
strain  at  yield. 


in  the  plastic  region  and  and  are  stress  and 


By  taking  o  and  e  at  4%  true  strain  and  assuming  the  stress- 
strain  curve  between  yield  and  4%  strain  is  linear  on  a  log-log 
basis,  the  following  strain  hardening  exponents  were  obtained: 
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Material 

Strain  Rate 
(Per  Sec) 

N 

(§  4%  Strain) 

Cu,  1/2-hard 

0.005 

0.038 

Cu,  1/2-hard 

950 

0.031 

Cu,  Ani.ealed 

0.005 

0.29 

Cu,  Annealed 

950 

0.27 

Strain  rate  behavior  is  summarized  in  the  cross  plots  given  in 
Figures  4  and  5.  The  as-received  material  shows  a  slight  rate 

3 

sensitivity  up  to  10  /sec  at  20°C,  while  the  yield  and  flow 
stresses  for  the  annealed  material  are  influenced  very  little 
(on  an  absolute  scale)  by  strain  irate.  At  elevated  temperatures, 
there  is  a  noticeable  increase  in  the  sensitivity  at  rates  above 
10/sec. 

Analysis  of  stran  rate  effects  under  uniaxial  stress  conditions 

(7) 

in  copper  by  Kumar  and  Humble  has  indicated  that  at  low  rates 
(<  10/sec)  flow  stress  is  controlled  by  thermal  activation  of 
mobile  dislocations,  primarily  over  forest  dislocation  barriers. 
At  higher  rates  (>  10^/sec),  dislocation  motion  becomes  viscous 
phonon  drag  limited.  They  also  concluded  that  mobile  disloca¬ 
tion  density  at  high  rates  is  strain  rate  independent. 

Although  deformation  twinning  in  shocJc- loaded  copper  (see  Equa¬ 
tion  of  State  discussion  below)  has  been  observed,  no  experi¬ 
mental  evidence  of  twinning  in  the  uniaxial  stress  tests  was 
(8) 

found.  Also,  plain-strain  compression  tests  on  annealed, 

(9) 

single-crystal  copper  by  Wonsiewicz  and  Chin  '  revealed  no 
evidence  of  twinning.  Their  tests  showed  large  strength  dif¬ 
ferences  (factor  of  three)  in  true  stress-true  strain  curves 
for  various  crystal  orientations.  However,  when  presented  as 
shear  stress-shear  strain  curves  and  corrected  for  crystal  ro¬ 
tation,  the  curves  overlapped  in  a  relatively  narrow  band. 
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Figure  5  Compressive  Yield  Stress  vs  Log 
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All  data  presented  in  Figures  3  to  5  are  for  a  test  direction 
parallel  to  the  bar  axis.  A  few  tests  conducted  normal  to  the 
axis  showed  the  flow  stress  to  be  'v-  0.25  kbar  higher,  with  rate 
sensitivity  and  strain  hardening  behavior  about  the  same. 


REVERSE  LOADING  TESTS 


it 

The  reverse  loading  behavior  or  Bauschinger  effect  was  studied 
by  performing  uniaxial  tensile  stress  tests  after  the  material 
had  been  prestrained  in  the  opposite  (compressive)  sense.  Two 
variations  of  this  method  of  determining  the  Bauschinger  effect 
were  used: 

1.  Uniaxial  stress  specimens  were  loaded  in  compres¬ 
sion  to  strains  of  0.5  to  3.0%,  unloaded  to  zero 
stress,  and  then  immediately  loaded  in  tension. 

2.  Specimens  were  compression- tension  cycled 
several  times  with  a  total  plastic  strain  of 
1%  per  cycle. 

The  materials  were  deformed  at  room  temperature  and  at  a  strain 

-3 

rate  of  10  /sec  in  the  above  experiments.  The  results  of  the 
first  tests  are  shown  in  Figure  6.  The  Bauschinger  strain 
(defined  here  as  the  plastic  strain  in  the  reverse  loading  path 
at  3/4  of  the  initial  yield  strength)  increases  sligh^^''”  ’  in¬ 
creasing  prestrain.  No  evidence  of  any  plastic  deformation  dur¬ 
ing  unloading  was  seen.  However,  plastic  deformation  during 
unloading  has  been  observed  at  strains  less  than  0.5%  with  mea¬ 
surement  sensitivities  much  greater  than  those  in  the  present 
experiments . 


The  Bauschinger  effect  is  the  reduction  in  stress  at  a  given 
strain  in  the  reverse  direction  compared  to  that  in  the  orig¬ 
inal  loading  direction. 
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The  effect  of  additional  strain  limited  cycles  is  shown  in 
Figure  7  with  the  material  work-softening  as  the  number  of 
previous  loading  cycles  increases  and  approaching  some  equili- 
britun  curve.  This  is  similar  to  that  seen  by  Morrow who 
found  many  materials  tend  to  approach  some  saturated  curve 
either  by  work-hardening  or  work- softening. 

Although  all  of  the  present  tests  were  performed  at  low  strain 
rates  (tests  in  the  order  of  seconds)  compared  to  those  seen 
during  shock  wave  loading,  the  data  should  give  some  indication 
of  material  behavior  during  reverse  loading.  The  cross-effect 
(influence  of  unloading  in  some  direction  other  than  exactly 
opposite  the  loading  path)  was  found  important  by  Hecker^^^^ 
and  others,  but  the  present  work  describes  the  behavior  on  only 
one  axis  on  the  biaxial  stress  plane.  For  accuracy  in  material 
modeling,  these  effects  should  be  considered  if  the  material 
model  must  be  used  for  a  complex  stress  deformation. 

Relating  the  reverse  loading  behavior  to  a  dislocation  mechanism 
is  difficult  in  the  uniaxial  stress  case  and  becomes  more  diffi¬ 
cult  as  the  stress  state  becomes  more  complex.  Trying  to  relate 
the  absolute  magnitudes  or  Bauschinger  strain  to  the  microstruc¬ 
ture  of  the  material  is  impossible  at  this  time  and  extrapolat¬ 
ing  the  present  data  to  materials  having  a  microstructure  dif¬ 
ferent  than  the  1/2-hard  OFHC  copper  tested  would  be  unjustified. 
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SECTION  III 

ELASTIC  BEHAVIOR 


Measurement  of  elastic  constants  provides  additional  inputs  to 
the  study  of  material  response.  The  three  basic  measurements 
are  density,  longitudinal  (dilatational)  wave  velocity  and 
shear  (transverse)  wave  velocity.  Values  obtained  at  standard 
temperature  and  pressure  (20°C  and  atmospheric  pressure)  estab¬ 
lish  the  parameters  used  in  first-order  (temperature  and  pres¬ 
sure  independent)  calculations  of  material  response.  Measure¬ 
ments  of  temperature  and  pressure  dependence  permit  more  exact 
calculations  to  be  made  and,  in  particular,  lead  to  prediction 
of  pressure-compression  isotherms,  adiabats  and  hydrostats.  If 
these  measurements  are  made  with  sufficient  accuracy  (<  0.1%), 
the  calculated  pressure-compression  behavior  can  frequently 
be  extrapolated  with  reasonable  confidence  to  higher  pressures 
than  covered  in  the  actual  measurements. 

WAVE  VELOCITY  MEASUREMENTS 

Initial  measurements  of  longitudinal  and  shear  wave  velocities 
were  made  at  20’C  and  atmospheric  pressure  (i.e.,  P  *  0) . 

Values  obtained  were: 

C_  »  4,757  ±  0.004  mm/usec 

ii 

Cg  »  2.247  ±  0.002  mm/usec 
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Measurements  were  then  made  of  the  tenperature  and  hydrostatic 
pressure  dependence  of  the  wave  velocities.  Correction  of  the 
data  for  thermal  expansion  and  hydrostatic  compression  gave  the 
following: 

Temperature  Dependence,  0-100®C,  P  =  0 
=  4.766  -  (0.000434  X  T)  mm/ysec 

‘Li 

C_  =  2.253  -  (0.000300  X  T)  mm/ysec 
o 

(T  in  “O 

Pressure  Dependence,  0-9  kbar,  T  =  20 ®C 
=  4.757  +  (0.00685  X  P)  mm/ysec 

Ij 

Cg  =  2.247  +  (0.00243  X  P)  mm/ysec 
(P  in  kbar) 

Within  the  limits  of  experimental  accuracy,  the  second-order 
coefficients  were  zero. 


Third-order  elastic  constants  can  be  determined  by  measuring 
elastic  wave  velocity  as  a  function  of  uniaxial  stress.  A 
copper  specimen  was  subjected  to  a  static  compressive  load  and 
C.  was  measured  normal  to  the  loading  direction.  Over  the 
stress  range  of  0  to  2  kbar,  no  change  in  Cj^  was  observed 
(after  correcting  for  lateral  strain)  within  the  accuracy  of 
the  measurements  {i  0.1%).  The  third-order  elastic  constants 
for  copper  single  crystals  have  been  reported  by  Salama  and 
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ELASTIC  CONSTANTS 


Assuming  the  material  is  isotropic  (C^  and  Cg  are  independent 
of  direction  of  measurement  in  a  polycrystalline  specimen) , 
wave  velocity  data  are  used  to  calculate  various  elastic  con¬ 
stants  At  20®C  and  zero  pressure,  the  following  adia¬ 

batic  constants  were  obtained  (see  Appendix  A) : 


Longitudinal  wave  velocity,  Cj^  =  4.757  mm/usec 

Shear  wave  velocity,  Cg  =  2.247  mm/usec 

Bulk  wave  velocity,  C^  =  3.99  iran/usec 

Sound  wave  velocity,  Cg  =  3.70  mm/usec 

Rayleigh  wave  velocity,  =  2.11  mm/usec 

Poisson's  ratio,  u  =  0.356 

Bulk  modulus,  K  »  1418  kbar 

Shear  modulus,  G  =  450  kbar 

Elastic  modulus,  E  »  1221  kbar 

Lame's  parameter,  A  >  m?  kbar 


The  bulk  and  shear  moduli!  are  used  in  the  development  of  the 
ultrasonic  equation  of  stat.^.  They  are  obtained  first  as 
functions  of  T  and  P  from  the  following: 


^  r  2 


A  cubic  material  is  elastically  isotropic  when 


'*11  '*12 

where  and  are  elastic  stiffness  constants. 

Although  single  crystal  copper  is  highly  anisotropic  with 

? — ^ —  *  s 764-n'26  “  this  anisotropy  may  be 


^ll“^12 


averaged  out  if  the  grains  are  randomly  oriented  crystal- 
lographically  and  the  strain  is  unifoi'm. 
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This  gave: 

Teinperature  dependence,  P  »  0 

Kg  =  1424  -  (0.278  X  T)  kbar 
Gg  =  453  -  (0.144  X  T)  kbar 

Pressure  dependence,  T  =  20“C 

Kg  =  1418  +  (5.53  :<  P)  kbar 
Gg  ^  450  +  (1.30  X  P)  kbar 

where  the  superscript  S  indicates  adiabatic. 

Finally,  isothermal  values  of  K  and  G  as  well  as  adiabatic 
and  isothermal  pressure  and  temperature  derivatives  at  20“C 
and  zero  pressure  were  calculated  using  the  equations  given 
in  Appendix  A.  Results  are  given  in  Table  I.  Use  of  these 
constants  in  calculating  the  isotherm,  adiabat  ard  hydrostat 
is  discussed  in  the  section  on  Ultrasonic  Equation  of  State. 


GRUNEISEN  PARAMETER ' ^ ^ ^ 

The  Cruneisen  ratio  y  is  a  parameter  originally  given  by 
Gruneiscn  in  a  solid  equation  of  state  relating  pressure  to 
volume  and  energy  (y  “  V  (?P/3E)y) .  This  parameter  can  be 
defined  thermodynamically  as: 
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TABLE  1 

ELASTIC  CONSTANTS  FOR  COPPER 
(P  =  0,  T  =  2J“C) 


PARAMETER 

VALUE 

k" 

1418  )cbar 

‘ll"^T  *  •'OT 

5.53 

-0.278  )cbar/“C 

,3K^.  „S' 

‘ir->s  =*  ^OS 

5.40 

,3K® 

<5t"’s 

12.1  kbar/^C 

1373  kbar 

,3k''^.  „T’ 

'3P  'T  “  "OT 

5.59 

3K^ 

'?I~>P 

-0.408  kbar/*C 

,3K^.  „T' 

^??">s  '  •'os 

5.42 

(1^-) 

12.1  kb«r/*C 

-  r/ 

450  kbar 

<If>t  -  Sr 

1.3C 

-0.144  kbar/*C 

“  Ss 

i.:4 

#s 

2,76  kbar/*C 
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For  an  isotropic  elastic  solid,  this  can  be  expressed  as 

^  ^  -  I  s')  '=> 

wh3re  Ct  and  C-  are  measured  under  adiabatic 

li  D 

conditions . 

At  20°C  and  zero  pressure  this  gives 

Yo=2.i3 

The  zero-pressure  Gruneisen  parameter  can  also  be  estimated 

(18) 

by  several  other  nethods,  including  those  of  Slater  , 

(19)  (20  21) 

Dugdale  and  MacDonald  ,  Anderson  and  Dienes  '  and 

(20  22) 

Schuele  and  Smith  '  .  These  methods  are  discussed  in 

Appendix  B. 

An  estimate  of  the  temperature  dependence  of  the  Gruneisen 
parameter  can  be  obtained  by  differentiating  Equation  4 
with  respect  to  temperature.  At  constant  pressure  (P  =  0, 

T  =  20°C) : 


The  pressure  dependence  can  be  estimated  by  differentiating 
with  respect  to  pressure  (see  Appendix  C) .  At  constant 
temperature  (T  =  20°C,  P  =  0) ; 
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Another  approach  is  to  assume  y/V  constant,  which  permits 
Oy/oP)^  to  be  estimated  directly  from  the  bulk  modulus 
data  (see  Appendix  C) .  At  constant  temperature  (20“C) : 


_ -7.76 _ 

(1373  +  5.59P)^*^^® 


=  -  0.0016Abar  @  P  =  0 
=  -  0.0010/kbar  @  P  =  100  kbar 
(23) 

McQueen  and  Marsh  used  shock  velocity  data  to  calculate 
Y  (assuming  y  =  f (V) ) ,  which  gave  Ay/AP  ~  -0.002  at  P  =  0  and 
Ay/AP  =  -0.001  at  P  =  100  kbar. 


DEBYE  TEMPERATURE 

The  Debye  temperature  6  is  a  parameter  related  to  the  maximum 
number  of  oscillation  frequencies  a  system  assumed  to  be  an 
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elastic  continuvun  can  support.  It  is  an  important  parameter 
in  thermal  energy  calculations  and  indicates  the  temperature 
above  which  variations  due  to  temperature  for  some  thermodynamic 
parameters  such  as  specific  heat,  thermal  expansion  and  Gruneisen 
ratio  become  small.  There  are  a  number  of  methods  for  calculat¬ 
ing  the  Debye  temperature,  and  one  approximation  suitable  for 
use  with  elastic  wave  velocity  data  for  polycrystalline  metals 
is : 


where  h  is  Planck's  constant,  k  is  Boltzmann's  constant, 
N  is  number  of  mass  points  and  V  is  sample  volume. 


At  20*C,  the  elastically-determined  Debye  temperature  for  copper 
is : 


O  <  •.  »  (tl\  1/3/1  .2  \-l/3 

6  =  4.3X10 

9  =  331^K 


(i^) 


for  N  =  4  (atoms  per  unit  cell),  V  =  (3.615  X  10  ^  mm) ^ 
(lattice  constant),  C  =  4.757  mm/ysec  and  C_  =  2.247  mm/ysec. 


This  compares  well  with  the  value  9  =  331. 3°K  obtained  by  a 
more  exact  computer  method  using  single-crystal  elastic 
data.'”' 
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SECTION  IV 

EQUATION  OF  STATE 


The  development  of  material  response  models  and  computational 
codes  requires  extensive  data  on  material  behavior  under  uni¬ 
axial  strain  conditions,  including  equation  of  state,  wave 
propagation  and  spall  fracture.  Most  data  are  used  directly 
in  developing  models  of  material  behavior,  although  independent 
check  data  (attenuated  wave  shapes  and  spall  wave  profiles)  are 
necessary  to  determine  accuracy  of  the  calculations.  This  sec¬ 
tion  presents  the  equation  of  state  analysis,  and  the  following 
two  sections  cover  wave  propagation  and  spall  fracture. 


SHOCK  WAVE  EQUATION  OF  STATE 

The  hugoniot  equation  of  state  is  the  locus  of  equilibrium 
states  reached  after  shocking  of  a  material,  and  the  data  are 
usually  obtained  either  as  stress-particle  velocity  points  or 
as  shock  velocity-particle  velocity  points.  The  shock  wave 
hugoniot  da«,_  for  copper  presented  in  this  report  were  obtained 
with  x-cut  quartz  gages  and  optical  techniques.  Representative 
quartz  gage  records  for  direct  impact  and  transmitted  wave  tests 
are  shown  in  Figure  8.  The  buffered  direct  impact  method  (tung¬ 
sten  carbide  buffer  plate  on  front  of  the  quartz)  permitted 
stresses  up  to  80  kbar  in  the  copper  while  keeping  the  stress 
in  the  quartz  at  an  acceptable  level.  The  use  of  quartz  gages 
for  transmitted  wave  tests  was  primarily  for  the  study  of  com¬ 
pressive  wave  development  and  elastic  precursor  decay  and  re¬ 
sults  are  discussed  in  the  Wave  Propagation  section.  The  direct 
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■DIRECT  IMPACT 

Cu  -►Q 

V  j-  0.  0906  mm/  /is 
9.  9  kbar 

■  DIRECT  IMPACT,  BUFFERED 

Cu-*-WC/Q 
V l“  0. 159  mm//is 
0",^=  44.0  kbar 

■  TRANSMITTED  WAVE 

AI-*-Cu/Q 
V i"  0. 165  mm/  /is 
18. 1  kbar 

Figure  8  Quartz  Gage  Records 
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impact  records,  which  give  a  stress-particle  velocity  point 
since  the  hugoniot  of  quartz  is  known,  showed  very  fast  rise 
times  {<  10  nsec)  and  then  a  constant  peak  stress  for  the  dura¬ 
tion  of  the  gage  transit  time.  The  slight  rounding  at  the  front 
of  the  direct  impact  records  in  Figure  8  is  due  to  impact  tilt, 
and  the  ramped  peak  stress  region  results  from  finite-strain 
effects  in  the  quartz  which  are  corrected  in  the  data  analysis. 

Hugoniot  data  at  stresses  of  100  to  200  kbar  were  obtained  by 
using  a  steak  camera  and  a  target  in  the  form  of  a  two-step 
cylinder  ("hat"  configuration) .  A  schematic  of  the  technique 
and  a  sample  record  are  shown  in  Figure  9.  The  two  surfaces 
of  the  target  are  viewed  through  the  camera  slit  and  shock  ar¬ 
rival  is  observed  as  a  change  in  surface  reflectance.  The 
average  shock  velocity  is  obtained  from  the  record,  and  par¬ 
ticle  velocity  is  assumed  to  be  1/2  Vj  (impact  velocity)  since 
impactor  and  target  are  the  same  material. 


The  hugoniot  may  be  expressed  in  several  forms,  A  convenient 
form  for  experimental  work  is  the  hugoniot  centered  about  zero 
stress  and  particle  velocity,  as  established  by  a  least-squares 
fit  to  the  data  in  the  stress-particle  velocity  (a„  -  u  )  plane. 

H  p 

Transformation  of  the  hugoniot  into  various  planes,  such  as 

shock  velocity-particle  velocity  (U  -  u  )  and  stress-volume 

S  P 

-  v) ,  is  performed  by  assuming  a  material  model.  This  was 
done  by  assuming  an  ideal  elastic-plastic  wave  structure  with 
equilibrium  initial  and  final  states  and  applying  the  mass  and 
momentiim  conservation  equations.  This  leads  to: 


a  + 
e 


(U  -  u^)  (u^  -  u^) 

S  e  p  e 


(10) 
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IMPACTOR 


TARGET 


CAMERA  SLIT 


Figure  9  Streak  Camera  Record,  "Hat"  Target 
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and 


(11) 


where  and  are  density  and  particle  velocity 

at  the  elastic  limit  a  . 

e 

Therefore,  when  either  a„  -  u  or  U_  -  u  relations  are  estab- 

H  p  S  p 

lished  and  a  ,  p  ,  u  and  are  known,  a  ,  U  ,  u  and  v  are 

6  6  0  Li  nop 

uniquely  determined.  The  resulting  equations  of  state  for 
copper  are  given  in  Figures  10,  11  and  12  and  are  listed  below: 


o„  =  1.5  +  344  u  +  143  u  ^  (Std.  Error=0.4  kbar) 
H  p  p 

a„  =  3.1  +  1308  u+3520 

n 

2 

Pjj  =  1308  y+3520  y  (Kydrostat) 

U-  =  3.915  +  1.53  u 
S  p 

where  y  =  (v^/v)-l 


(12) 


These  equations  of  state  are  based  on  the  measurements  indi¬ 
cated  by  the  symbols  in  Figure  10  and  =  8  kbar,  and  are 
considered  to  be  valid  over  the  stress  range  of  10  to  200  kbar. 
Accuracy  of  the  data  points  are  estimated  to  be  ^  2%  for 

and  .t  1%  for  u  . 

P 

The  hydros tat  (P  -y)  was  obtained  by  subtracting  the  stress- 

H 

axis  intercept  value  from  the  equation.  This  approach 
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satisfies  the  requirement  that  the  P„-)j  function  pass  through 

n 

the  origin,  assumes  che  deviatoric  stress  (a^)  is  constant 
(no  strain  hardening),  and  implies  a  value  of  3.1  kbcr  for 
Op.  Alternative  ways  of  estimating  Oj^  are  discussed  in 
Appendix  D. 

The  low-pressure  hugoniot  of  copper  (99.95+%)  has  been  reported 

iyfi\ 

previously  by  Munson  and  Barker, ^  '  who  used  the  flat-plate 

impact  technique  and  free  surface  velocity  measurements  to 
determine  the  hugoniot  up  to  116  kbar.  They  calculated  the 
hydrostat  (hugoniot  mean  pressure) ,  adiabat  and  isotherm,  using 
corrections  for  the  deviatoric  stress,  entropy  change  and  tem¬ 
perature  change.  The  hydrostat  was  given  as: 

=  0.734  X  10"^P„  -  1.41  X  10"S„^  (13) 

n  H 

Comparing  this  with  the  hydrostat  obtained  in  the  present 
study  gives: 


y 

P„  (Eq.  12) 
kbar 

P„  (Eq.  13) 
kbar 

0.020 

27.6 

28.2 

0.040 

58.0 

59.1 

0.060 

91.2 

94.1 

The  agreement  is  reasonable  below  100  kbars . 

The  high-pressure  hugoniot  of  copper  was  recently  reported 
by  Isbell,  et  as: 


=  3.964  +  1.463  u  (14) 

S  p 
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This  relation  is  based  on  shock  velocity-particle  velocity 
measurements  over  the  range  of  1  to  5  nibar. 


ULTRASONIC  EQUATIONS  OF  STATE 

The  equation  of  state  can  also  be  determined  from  ultrasonics 
measurements  (see,  e.g..  Ref.  28  to  31).  The  very  afifeirate 
measurement  of  elastic  wave  velocities  (and,  therefor^^ bulk 
modulus)  as  a  function  of  hydrostatic  pressure  makes ji® possible 
to  directly  estimate  shock  wave  compression  behavior  c®to  com¬ 
bine  the  ultrasonics  data  with  shock  wave  data  to  pre-JTct  very 
high  pressure  response.  Numerous  analytical  and  empirical  re¬ 
lations  have  been  developed  for  relating  pressure#  volume  and 

{32  33)  {34) 

bulk  modulus,  including  those  Birch,'  '  '  Murnaghan'  ' 

and  Keane .  The  use  of  these  relations  is  discussed  in 

Appendix  E  and  the  results  are  s^unmarized  below; 

Birch  EOS 


=  2060 


(...) "  -  (...) 


5/3 


( 

/  \  2/3 

1+1.19  { 

i 

\  /  ) 

(15) 


Murnaghan  EOS ,  Isotherm 


Pt  =  245 


.6  [  (i.m)  5-59  -ij 


(16) 


Murnaghan  EOS,  Adiabat 


Pg  =  262 


.6  Ul.,)  S-"”  -l] 


(17) 


34 


MANUFACTURING  DEVELOPMENT  •  GENERAL  MOTORS  CORPORATION 


MSL-70-23,  Vol.  V 


Keane  EOS, 


3.7 


Pg  =  559 


.3 


651.5  In 


(...) 


(18) 


Keane  EOS,  K^'  =  2.73 


®|||||^  1025  -1383  In  ^1+u^ 


(19) 


In  the  cLbove  equations,  subscripts  T  and  S  indicate  isothermal 
and  isentropic,  respectively. 


Duvall  gives  a  correction  for  the  change  in  entropy  across  a 
shock  front  based  on  the  Murnaghan  equation  of  state: 


^  (‘'os  *  l)  (l^) 


'■h  =  ^'s  ^  ' 


(20) 


Zel'dovich  gives  an  expression  relating  ®  and  P  ,  assuming 

/  .  ^  (30,38)  ^  ^ 

Y/v  is  constant:  ' 


Pg  =  Pt  t  3TqK 


o[-p|’^o(ifr)}-  H] 


Pg  =  P.J  +  20.8 


exp  {  2.13 


(w)  I -(.•.)] 


(21) 
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The  above  equations  of  state  are  con^ared  to  the  shock  wave 
EOS  in  Table  II.  The  entropy  correction  based  on  the  Murnaghan 
equation  has  been  applied  to  the  Birch  and  Keane  equations  as 
well.  The  ultrasonic  equations  of  state  show  good  agreement 
with  each  other  but  all  give  mean  pressure  values  at  a  given 
compression  that  are  4-8%  higher  than  the  shock  wave  hydros tat, 
up  to  250  kbars. 


The  shock  velocity-particle  velocity  relationship  can  also  be 
determined  from  ultrasonic  data.  Following  the  method  of 


Ruoff 


one  obtains: 


U„  =  C„  +  Su^  +  Au  ' 
S  B  p  p 


where 


■{ 


S  =  1/4 


(■'os  ^  l) 


For  copper: 


24  C, 


{’  -  4's 


+  4y, 


c  c  ” 


=  3.99  +  1.60  u  +  0.17  u 
S  P  P 

2  S " 

The  Up  coefficient  (0.17)  was  calculated  assuming  K^g  =  0. 

A  positive  second-order  coefficient  implies  a  concave  upward 

U„  -  u  curve.  However,  comparison  of  low-pressure  (Eq.  12) 
o  P 

and  high-pressure  (Eq.  14)  shock  wave  data  indicates  that  the 
slope  of  the  Ug  -  u  curve  decreases  slightly  with  increasing 
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h 

pressure.  If  one  assumes  that  the  relation  given  above  for 

S  " 

'A'  is  valid  then  Kgg  can  be  estimated  as  a  function  of  'A' . 

For  example,  for  a  linear  0_  -  u  relation,  A  =  0  and 
<?"  ®  P  /AA) 

=  -0.0057.  Pastine  and  Piacesi '  '  calculated  A  =  0.0144 

(23)  S" 

(using  data  of  McQueen  and  Marsh'  ’)  which  gives  =  -0.0062. 

s  ** 

It  would  seem  that  K^g  =  -0.006  for  copper  at  zero  pressure, 
which  gives  the  following  as  a  more  reasonable  ultrasonic 
Ug  -  Up  relation  for  low  pressures: 

=  3.99  +  1.60  u  -  0.01  u  ^  (24) 

S  p  p 

This  is  compared  with  the  linear  U„  -  u  relation  obtained 

S  p 

from  low-pressure  shock  wave  data  in  Figure  12. 


YIELD  BEHAVIOR 

A  complete  description  of  the  low-pressure  equation  of  state 
requires  consideration  of  yield  behavior  or  the  hugoniot  elas¬ 
tic  limit.  Values  of  compressive  yield  were  obtained  by  three 
independent  methods. 

A.  Uniaxial  Stress 

The  yield  level  was  determined  in  uniaxial  stress 
(0^)  as  a  function  of  strain  rate  and  then  con¬ 
verted  to  uniaxial  strain  using: 


This  comparison  assumes  no  phase  change  and  no  melting. 
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strain  Rate 
-1 

sec 


10' 

10 


-3 


e 

kbar 


7.0 

5.8 


B.  Wave  Profiles 

The  elastic  precursor  level  was  measured  as  a  function 
of  propagation  distance  using  quartz  gages .  These 
profiles  are  discussed  in  the  section  on  Wave  Propaga¬ 
tion. 


Propagation 

a 

e 

Distance 

mm 

kbar 

2 

'\/4.0 

12 

'\^2.8 

C.  Hugoniots 

The  yield  level  in  uniaxial  strain  was  inferred  by 

comparison  of  the  elastic  (PnC, )  and  plastic 

(Ojj  vs.  Up)  hugoniots,  i.e.,  by  defining  the  elastic 

limit  as  the  intersection  of  the  elastic  response 

line  and  the  curve-fit  to  the  stress-particle  velocity 

data: 

10  PrtC-U  =  A  +  Bu^  +  Cu^^ 

0  L  e  e  e 

solve  for  u 

e 

then  =  10  p-C  u 
e  u  L  e 
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Ug  =  0.019  min/jis 
=  8.0  koar 

G 

For  an  isotropic,  poly crystalline  metal  with  elastic-perfectly 
plastic  behavior  and  no  strain-rate  or  time-dependent  effects, 
one  would  expect  the  same  yield  level  to  be  obtained  by  each 

r 

I  method.  However,  copper  shows  both  strain- rate  sensitivity 

[  and  decay  of  the  elastic  precursor.  Also,  the  elastic  wave  is 

I 

^  not  well-defined  and  shows  spreading  (increasing  rise  time)  as 

well  as  decay  with  distance.  Because  of  these  effects,  the 
elastic  limit  is  not  single-valued  and  a  compromise  must  be 
made  in  selecting  a  value  for  use  in  the  hugoniot  transformation 
equations  given  earlier  where  elastic-perfectly  plastic  behavior 
was  assumed. 

The  elastic  limit  calculated  from  the  hugoniots  (8  kbar)  was 

used  in  determining  the  a„  -  y  and  U  -  u  equations  given  for 

"  S  P 

copper  in  the  low-pressure  region.  Extrapolation  of  the  uni¬ 
axial  stress  tests  to  10^/sec  would  give  an  equivalent  yield 
in  uniaxial  strain  of  8  kbar.  Although  this  extrapolation 
is  very  uncertain,  the  strain  rate  in  the  elastic  wave  front 
in  copper  is  about  30^  to  10®/sec.  Therefore,  the  use  of  8  kbar 
as  the  hugoniot  elastic  limit  is  consistent  with  the  uniaxial 
stress  measurements.  The  yield  as  determined  from  the  wave  pro¬ 
files,  however,  is  considerably  lower  than  8  kbar,  at  least  for 
propagation  distances  greater  than  2  mm.  Also,  the  precursor 
amplitude  attenuates  >  80%  in  copper  in  the  first  2  ram,  which  is 
much  higher  than  observed  for  other  metals .  The  reasons  for 
this  apparent  disagreement  are  not  known,  but  may  be  related  to 
some  microscale  metallurgical  factor  such  as  microyielding  or 
to  the  elastic  anisotropy  of  single  crystal  copper. 
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METALLOGRAPHIC  STUDIES 

Grain  structure  in  copper  subjected  to  shock  loading  is  shown 

in  Figure  13.  Compressive  pulse  duration  for  these  tests  was 

lysec.  There  is  some  evidence  of  macrotwinning  although  it 

has  affected  a  relatively  small  volume  and  did  not  increase 

appreciably  from  18  to  86  kbar,  and  the  increase  in  hardness 

(41) 

was  small  ('^<  78  to  84  Rp) .  Smith  found  the  macrotwin  den¬ 
sity  to  still  be  low  at  225  kbar  but  to  become  substantial  at 
425  kbar.  It  is  unlikely  that  the  elastic  wave  behavior  in 
copper  can  be  attributed  to  the  effects  of  macrotwin  formation, 
especially  when  there  is  evidence  that  the  twinning  may  occur 
during  the  release  portion  of  the  compressive  pulse. 

Shock-loaded  copper  also  shows  twinning  on  a  microscale,  i.e., 

twins  on  the  order  of  0.5  ym  width  which  are  too  small  to  be 

seen  optically.  Examples  from  transmission  electron  microscope 

of  shock-loaded  copper  foil  (not  in  the  same  initial  condition 

as  the  bar  stock)  are  shown  in  Figure  14.  The  dislocation  den- 

10  2 

sity  in  the  as-received  foil  was  estimated  at  10  /cm  with 
a  few  microtwins  observable,  while  the  shocked  (34  kbar, 

0.5  ysec  pulse)  foil  appeared  to  have  a  higher  dislocation 
density  with  definite  indication  of  increased  microtwinning. 
Transmission  electron  microscope  studies  by  other  investigators 
(see  e.g..  Refs.  42  to  46)  have  indicated  that  microtwinning 
may  occur  at  stresses  at  low  as  16  kbars,  and  becomes  substantial 
in  the  range  of  100  to  200  kbars.  These  studies  also  showed  sub¬ 
stantial  dislocation  cell  structure  development  with  dislocation 
densities  on  the  order  of  10  /cm"  for  stress  levels  of  several 
hundred  kbars. 
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Figure  13  Shock  Wave  Test  Specimens,  Copper 
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MEl ALLOGRAPHIC  STUDIES 

Grain  structure  in  copper  subjected  to  shock  loading  is  shown 

in  Figure  13.  Compressive  pulse  duration  for  these  tests  was 

'V'  lysec.  There  is  some  evidence  of  macrotwinning  although  it 

has  affected  a  relatively  small  volume  and  did  not  increase 

apprec'ably  from  18  to  86  kbaur,  and  the  increase  in  hardness 

(41) 

was  smalJ  78  to  84  Rp) .  Smith  found  the  macrotwin  den¬ 
sity  to  still  Le  low  at  225  kbar  but  to  become  substantial  at 
425  kbar.  It  is  unlikely  that  the  elastic  wave  behavior  in 
■opper  can  be  attributed  to  the  effects  of  macrotwin  formation, 

especially  when  there  is  evidence  that  the  twinning  may  occur 

(41  42) 

during  the  re  ^ase  portion  of  the  compressive  pulse.  ' 

Shock-loaded  copper  also  shows  twinning  on  a  microscale,  i.e., 

twins  on  the  order  of  0.5  pm  width  which  are  too  small  to  be 

seen  optically.  Examples  from  transmission  electron  microscope 

of  shock-loaded  copper  foil  (not  in  the  same  initial  condition 

as  the  bar  stock)  are  shown  in  Figure  14.  The  dislocation  den- 

10  2 

sity  in  the  as-received  foil  was  estimated  at  ■v  10  /cm  with 
a  few  microtwins  observable,  while  the  shocked  (34  kbar, 

0.5  usee  pulse,  foil  appeared  to  have  a  higher  dislocation 
density  with  definite  indication  of  increased  microtwinning. 
Transmission  ele'^tron  microscope  studies  by  other  investigators 
(see  e.g..  Refs.  42  to  46)  have  indicated  that  microtwinning 
may  occur  at  stresses  at  low  as  16  kbars,  and  becomes  substantial 
in  the  range  of  100  to  200  kbars.  These  studies  also  showed  sub¬ 
stantial  dislocation  cell  structure  development  with  dislocation 

11  2 

aensitxes  on  the  order  of  10  /cm  for  stress  levels  of  several 
hundred  kbars. 
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SECTION  V 

WAVE  PROPAGATION 


Compressive  wave  development  and  elastic  precursor  decay  are 
used  as  basic  input  data  for  material  response  models,  while 
a  complete  understanding  of  shock  wave  propagation  requires 
inclusion  of  release  wave  behavior.  Wave  attenuation  data 
provide  the  most  stringent  check  on  the  accuracy  of  shock 
wave  response  calculations . 


COMPRESSIVE  WAVE  BEHAVIOR 

Profiles  of  shock  waves  propagated  through  a  specimen  are  ob¬ 
tained  as  stress-time  or  velocity-time  histories.  Profiles  may 
be  used  in  this  form  directly  for  comparison  with  calculated  or 
predicted  material  response.  When  it  is  necessary  to  transform 
the  data  to  obtain  material  stress  from  transducer  stress  or 
particle  velocity  from  interface  velocity,  an  impedance  matching 
technique  is  applied,  utilizing  the  conservation  relations  and 
assuming  time-independent  behavior.  The  implications  of  this 
assumption  must  be  considered  when  utilizing  the  transformed 
data. 

Structure  in  the  compressive  wave  front  in  copper  is  shown  in 

★ 

the  quartz  gage  data  in  Figure  15.  The  wave  front  is  character¬ 
ized  by  a  non- steady- state  elastic  portion  and  a  transition  to 
a  spreading  plastic  wave  (nominal  final  stress  of  18  kbar) .  The 
elastic  wave  would  be  overdriven  by  the  plastic  wave  at  200  kbar. 

* 

In  this  and  subsequent  figures  containing  quartz  gage  data, 
tilt  refers  to  the  time  required  for  a  step-input  to  sweep 
across  the  gage  electrode  diameter. 
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Figure  15  Compressive  Waves,  Propagation  Distance 
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The  rise-time  of  the  elastic  wave  is  too  large  to  be  attributed 
to  tilt  effects.  It  seems  evident  that  only  the  initial  elastic 
disturbance  is  traveling  with  longitudinal  wave  velocity  (C  ) , 
and  subsequent  stress  increments  in  the  elastic  wave  are  travel¬ 
ing  with  lower  velocity.  This  behavior  is  not  appreciably  in¬ 
fluenced  by  peak  stress,  as  shown  in  Figure  16.  The  front  por¬ 
tion  of  the  waves  overlap  almost  exactly  for  peak  sti esses  of 
18  and  44  kbar.  Compressive  waves  in  1/2-hard  and  annealed 
copper  are  compared  in  Figure  17.  The  elastic  portion  has  almost 
disappeared  for  the  annealed  material  Jjs&e  Figure  2  for  micro- 
structure) .  These  profiles  indic^jb^  that  wave  structure  is  very 
sensitive  to  thermal-mechanica^T^istory  of  the  specimen. 


Although  the  conditj^ohs  for  development  of  steady-state  plastic 
waves  in  copper^^were  not  studied  in  detail  (i.e.,  the  minimum 
propagation  distance  for  a  given  stress  level  at  which  the 
pleas tic  wave  rise- time  becomes  constant  was  not  established) , 
the  rise-time  in  the  plastic  wave  increases  as  propagation  dis¬ 
tance  goes  from  5  to  10  mm  at  18  kbar  (see  Figure  15) .  Also 
shown  in  Figures  15  and  16  are  the  shock  wave  velocities  as 
calculated  from  the  shock  wave  hugoniot  (Equation  12) .  Al- 
through  there  is  qualitative  agreement  in  that  the  predicted 
plastic  wave  velocity  lies  within  the  rise-time  of  the  measured 
wave,  the  plastic  wave  cannot  be  considered  to  have  a  unique 
velocity  for  a  given  stress  level,  at  least  for  stresses  below 
50  kbar.  This  means  that  transformation  of  the  hugoniot  from 

the  a„  -  VI  plane  to  a„  -  u  or  U_  -  u  by  assuming  ideal  elas- 
n  p  n  op 

tic-plastic  v;ave  structure  and  steady-state  conditions  can  lead 
to  considerable  uncertainty  at  lower  stresses. 

Elastic  precursor  structure  in  copper  has  also  been  studied  by 
(47) 

Jones  and  Mote  using  flat-plate  impact  techniques  and 
quartz  gages.  For  annealed  single  crystals,  they  found  elastic 
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waves  of  1  to  2  kbars  amplitude  (depending  on  crystal  orienta¬ 
tion)  ,  with  fast  rise  time  ('v-  10  nsec),  and  stress  relaxation 
after  the  elastic  wave.  However/  for  prestrained  (3-1/2%) 
single  crystals  as  well  as  polycrystalline  copper  they  reported 
essentially  zero  yield  at  the  wave  front,  with  the  "elastic 
wave"  being  a  ramped  transition  to  the  plastic  wave.  This 

behavior  was  attributed  to  differences  in  dislocation  densities, 

6  2  9  2 

being  estimated  at  10  /cm  for  the  annealed  material  and  10  /cm 

for  the  prestrained  material. 

Although  the  other  PREDIX  metals  (6061-T6  aluminum,  alpha  titan¬ 
ium  and  tantalum)  showed  elastic  wave  risetimes  less  than  10  nsec, 
a  ramped  elastic  wave  has  been  observed  in  polycrystalline  beryl- 
lium  and  uranium.  '  Although  polycry-'  talline  beryllium  has 
a  dispersive  elastic  wave,  single-crystal  beryllium  shows  a  fast¬ 
rising  elastic  wave  with  about  a  factor  of  10  variation  in  pre¬ 
cursor  level  depending  on  crystal  orientation  (with  propagation 
along  the  C-axis  giving  the  highest  value) .  This  anisotropy 
could  contribute  to  smearing  of  the  wave  front  in  a  polycrys tal¬ 
line  material  with  randomly  oriented  grains. 

Dislocation  models  for  predicting  compressive  wave  development 
utilize  measurements  of  elastic  precursor  decay  to  evaluate 
various  constants  related  to  dislocation  density  and  multipli- 
cauxon. Unfortunately,  it  is  very  difficult  to  obtain 
accurate  elastic  precursor  level  data  at  propagation  distances 
less  than  2  mm  and  it  is  in  this  region  (0  to  2  mm)  that  most 
of  the  precursor  decay  occurs  for  many  metals  of  interest. 

Also,  even  at  larger  distances  there  is  always  some  uncer¬ 
tainty  in  the  precursor  level  due  to  both  material  scatter 
and  dispersion  characteristics  and  to  experimental  errors 
(particularly  tilt) .  The  result  is  that  much  of  the  data 
presently  available,  when  given  realistic  uncertainty  bars. 
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can  be  used  in  support  of  a  number  of  different  models  of 
precursor  attenuation.  Copper  is  no  exception  to  this,  as 
shown  in  Figure  18.  The  calculated  initial  elastic  impact 
stress  was  21  kbar  for  these  tests  and  the  precursor  level 
has  dropped  80%  to  4  kbar  after  only  2  mm  of  travel. 

The  data  for  X  >  2  mm  were  obtained  from  quartz  gage  records 
such  as  those  in  Figure  15.  Although  profiles  were  obtained 
for  propagation  distances  beyond  13  mm,  it  was  not  possible  to 
determine  a  precursor  level  at  these  distances  because  of  the 
excessive  dispersion  of  the  wave  front.  An  area-monitoring 
gage  such  as  quartz  is  obviously  more  sensitive  to  tilt  than 
a  point-monitoring  gage  such  as  the  velocity  interferometer 
which  shows  the  deviation  from  elastic  behavior  more  clearly 
(although  the  point-monitoring  gage  is  more  susceptible  to  an¬ 
isotropic  effects  since  the  spot-size  for  the  interferometer 
may  be  on  the  order  of  the  grain  size  in  the  material) .  The 
interferometer-measured  elastic  waves  did,  however,  show  rise- 
times  >  25  nsec  for  propagation  distances  >  2  mm.  One  velocity 
interferometer  test  with  a  0.52  mm  specimen  gave  a  precursor 
level  of  6  ±  1  kbar,  which  indicates  that  the  decay  curve  does 
start  to  break  upward  sharply  at  1-2  mm. 

The  uncertainty  bars  in  Figure  18  reflect  the  dispersion  and 
rounding  at  the  elastic  front.  This  spreading  is  due  to  at 
least  two  factors.  First,  since  the  quarts:  gage  averages  stress 
over  the  electrode  area,  small  differences  in  wave  front  arrival 
times  at  the  specimen/gage  interface  would  give  a  finite  rise¬ 
time  rather  than  an  instantaneous  stress  jump.  These  small 
differences  in  arrival  times  may  be  due  to  wave  reflections  and 
dispersion  at  grain  boundaries  and  to  slight  differences  in  local 
wave  velocity  arising  from  grain  anisotropy.  Second,  the  influ¬ 
ence  of  shock  wave  tilt  on  a  finite-area  gage  is  to  smooth  out 
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Figure  18  Elastic  Precursor  Decay  in  Copper 
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abrupt  changes  in  stress  level  as  well  as  to  increase  recording 
rise-time  in  the  wave  front.  This  becomes  particularly  severe 
at  short  propagation  distances  where  separation  of  the  elastic 
and  plastic  waves  is  only  on  the  order  of  tens  of  nanoseconds. 
The  effect  of  tilt  is  demonstrated  in  Figure  19.  For  a  very  low 
amount  of  tilt  (4  nsec  as  measured  across  the  quartz  gage  elec¬ 
trode)  there  is  a  distinct  transition  from  the  elastic  to  the 
plastic  wave,  while  for  tilt  of  31  nsec  the  structure  at  the 
front  of  the  wave  is  almost  completely  lost.  The  data  have  been 
plotted  to  give  the  same  breakaway  point  (corresponding  to  wave 
arrival  at  the  edge  of  the  electrode  area) ,  which  accounts  for 
the  offset  between  the  waves  in  the  plastic  region. 


RELEASE  WAVES  AND  WAVE  ATTENUATION 


Release  waves  ir.  copper  were  studied  by  using  relatively  thin 
copper  in^actors  and  measuring  the  conplete  wave  profile  with 
the  laser  velocity  interferometer.  After  impact,  elastic  and 
plastic  waves  propagate  from  the  impact  surface  into  both  im- 
pactor  and  target.  The  waves  that  reach  the  impactor  free  rear 
surface  reflect  as  rarefaction  waves,  unloading  the  material. 
These  unloading  waves  propagate  across  the  impact  surface  into 
the  target  and  are  recorded  at  the  target  rear  surface  after 
arrival  of  the  con^ressive  wave.  Tests  conditions  for  all 
velocity  interferometer  data  discussed  in  this  report  are 
listed  in  Table  III.  The  test  numbers  correspond  to  those 
listed  next  to  each  wave  profile  given  in  subsequent  figures. 

Unattenuated  wave  profiles  for  copper  are  shown  in  Figure  20. 

A  fused  quartz  window  was  used  on  the  rear  of  the  copper  target 
to  reduce  the  magnitude  of  the  unloading  due  to  reflection  of 
the  congress! ve  wave.  This  suppresses  spall  fracture  in  the 
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TABLE  III 


VELOCITY  INTERFEROMETER  TEST  DATA  (20®C) 


TEST 

NO. 

IMPACT 

VELOCITY 

(mm/ps) 

^0 

IMPACTOR 

THICKNESS 

(nun) 

X 

TARGET 

THICKNESS 

(nun) 

X/Xg 

CONFIG . 

149 

0.134 

0.930 

.  2.484 

2.67 

CU-CU 

150 

0.1136 

1.008 

2.502 

2.48 

Cu-Cu 

151 

0.217 

0.881 

2.499 

2.84 

Cu-Cu 

157 

0.318 

0.856 

2.413 

2.82 

Cu-Cu/FQ 

158 

0.217 

0.983 

2.512 

2.56 

Cu-Cu/FQ 

160 

0.228 

0.993 

5.075 

5.11 

Cu-Cu/FQ 

161 

0,227 

0.919 

10.117 

11.01 

Cu-Cu/FQ 

174 

0.244 

0.495 

3.785 

■ 

Cu-Cu 

0. 3i8mm/A»sec 


Figure  20  Wave  Profiles,  Stress  Dependence 
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target  which  would  prevent  observation  of  the  coirpressive 
release  wave  from  the  impactor.  The  measured  velocity  then 
is  interface  velocity,  which  is  greater  than  particle  velocity 
but  less  than  impact  velocity  since  the  inpedance  of  copper  is 
greater  than  fused  quartz.  Note  the  ranged  elastic  wave,  the 
rouiided  top  of  the  plastic  wave,  and  the  structxire  in  the  re¬ 
lease  wave.  This  release  wave  structure  is  due  to  two  effects. 
First,  elastic  and  plastic  waves  have  spread  in  the  impactor 
before  reaching  the  free  rear  surface  of  the  impactor  such 
that  all  the  unloading  does  not  occur  at  the  same  time,  i.e., 
the  unloading  or  release  wave  does  not  start  out  as  an  abrupt 
pressure  drop.  Secondly,  additional  structure  may  be  developed 
due  to  elastic-plastic  release  in  the  material.  It  is  difficult 
to  separate  these  two  contributions  to  release  wave  structure, 
however,  the  measured  profiles  do  serve  as  check  data  for  shock 
wave  calculations. 

If  target  thickness  is  great  enough  coirpared  to  impactor  thick¬ 
ness,  the  release  wave  will  overtake  the  compressive  wave  and 
attenuate  the  peak  pressure.  Wave  attenuation  is  shown  for 
copper  in  Figure  21  for  an  initial  stress  of  41  kbar.  The  over¬ 
taking  point  is  between  X/X^  =  5  and  10,  with  the  attenuation 
being  12%  at  X/Xq  =  11.  This  is  a  relatively  low  attenuation 
rate  compared  with  6061-T6  aluminum,  alpha  titanium  and  tantalum, 
which  show  35%,  35%  and  50%  attenuation  respectively,  at  X/X^ 

=  11  (for  comparable  stresses  and  thicknesses)  . 
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SECTION  VI 

SPALL  FRACTURE 


Spall  fracture  by  plate  impact  results  from  reflection  of 
compressive  waves  from  a  relatively  low  impedance  interface 
(normally  a  free  surface)  and  their  subsequent  interaction. 
Spall  studies  were  carried  out  with  both  passive  and  active 
techniques.  Passive  methods  involve  the  recovery  and  examina¬ 
tion  of  shock-loaded  specimens.  Metallographic  exeunination 
establishes  the  type  and  degree  of  damage,  which  can  be  cor¬ 
related  with  impact  parameters  such  as  velocity,  impactor 
thickness  and  target  thickness.  Active  methods  utilize  the 
laser  velocity  interferometer  and  provide  quantitative,  time- 
resolved  data  on  the  influence  of  internal  fractures  or  spall 
surfaces  on  shock  wave  profiles,  as  measured  at  the  rear  sur¬ 
face. 


RECOVERY  TESTS 


The  spall  behavior  of  copper  was  studied  by  carrying  out  a 
series  of  impact  and  recovery  tests,  where  the  target  was 
sectioned  across  a  diameter,  polished,  etched  and  examined 
optically  at  a  magnification  of  'v  SOX.  The  specimens  were 
then  graded  or  classified  according  to  the  degree  of  fracture 
that  was  observed,  which  ranged  from  no  visible  damage  to  com¬ 
plete  material  separation.  The  incipient  spall  threshold  vv’as 
defined  as  the  impact  velocity  (for  a  given  set  of  impact  par¬ 
ameters)  corresponding  to  cracking  over  at  least  50%  of  the 
width  of  the  section  estimated  to  be  under  a  condition  of  plane 
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strain  during  the  time  of  loading.  The  specimens  were  polished 
and  etched  to  provide  maximum  contrast  between  cracks  and  sound 
material.  The  procedure  followed  was: 

1.  Wet  grind,  320/400/600  grit  silicon  carbide  paper. 

2.  Rough  polish,  6  um  diamond  paste  on  rayon  cloth. 

3.  Final  polish,  1  ym  diamond  paste  on  rayon  cloth. 

4.  Etch,  equal  parts  water  (H2O) ,  hydrogen  peroxide 
(3%  H2O2)  and  ammonium  hydroxide  (NH^OH) , 

10-20  sec  swab  etch,  20‘*C. 

Optical  photomicrographs  of  recovered  specimens  are  shown  in 
Figures  22  and  23.  Fracture  is  of  a  ductile  nature  with  the 
development  of  individual  voids.  As  impact  velocity  increases, 
the  number  of  voids  becomes  greater  but  their  size  increases 
only  slightly.  At  high  enough  velocities,  the  voids  begin  to 
coalesce  and  a  complete  fracture  surface  is  formed.  Note  that 
as  velocity  increases,  the  voids  tend  to  cluster  within  a  very 
narrow  band  with  only  a  few  small  voids  outside  this  band.  Al¬ 
though  no  correlation  has  been  proven,  copper  has  a  poorly-de¬ 
fined,  low-amplitude  elastic  wave  and  shows  a  definite  "spall 
plane",  while  the  other  PREDIX  metals which  have  well- 
defined  elastic  wave  structures  show  considerable  dispersion 
of  the  spall  plane,  as  demonstrated  in  Figure  24.  This  behavior 
can  probably  be  related  to  differences  in  tensile  stress  history 
in  the  region  of  the  "nominal"  spall  plane.  For  copper,  the 

it 

average  distance  from  the  spall  plane  to  the  target  rear  surface 
was  approximately  equal  to  the  impactor  thickness,  'or  unatten¬ 
uated  pulses.  A  comparison  of  spall  fracture  in  sectioned  speci¬ 
mens  as  observed  by  optical  and  scanning  electron  microscopy  is 

* 

All  spall  photomicrographs  in  this  report  have  the  target 
rear  surface  at  the  top  of  the  picture. 
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Figure  24  Spall  Fractures  in  Metals 
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given  in  Figure  25.  Although  the  greater  depth  of  field  of 
the  SEM  gives  more  detail  in  and  around  the  voids  (particularly 
at  higher  magnifications)  the  optical  microscope  is  adequate 
for  general  examination  and  determination  of  the  spall  threshold. 

Scanning  microfractographs  of  the  fracture  surface  for  a  com¬ 
pletely  separated  specimen  are  shown  in  Figure  26.  The  equi- 
axed  nature  of  the  dimples  indicates  failure  by  a  normal  rup¬ 
ture  process/  where  the  voids  have  grown  and  coalesced  under  the 
action  of  plastic  strain  in  the  direction  of  wave  propagation. 
The  dinples  are  relatively  deep  and  the  walls  show  a  rippled 
structure  due  to  deformation  or  stretching.  These  ripples  may 
be  serpentine  glide  patterns  resulting  from  intersection  of  slip 
planes  with  the  void  w&lls  during  growth. 

Spall  data  for  all  tests  with  unattenuated  pulses  (X/X^  =  2) 
are  given  in  Figure  27.  The  symbols  are  defined  as: 

•  Complete  Separation 
9  Above  Incipient 

•  Incipient  Spall 
9  Below  Incipient 

O  No  Visible  Damage 

As  is  typically  the  case  with  metals,  the  impact  velocity  for 
incipient  spall  increases  for  decreasing  impactor  thickness. 

This  implies  that  the  peak  stress  (and,  therefore,  strain)  re¬ 
quired  to  create  spall  fracture  increases  with  decreasing  pulse 
width  or  time  of  loading.  The  incipient  spall  results  from  all 
tests,  including  an  attenuated  wave  series,  are  summarized  in 
Table  IV. 
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OPTICAL  MICROSCOPE  SCANNING  MICROSCOPE 


Figure  25  Spall  Fractures  in  Copper 

(1.5  nun  -*■  3.0  nun,  0.115  nun/ps) 


Figure  26  Spall  Fracture  Surface  in  Copper 
(0.5  Rin  i.O  BUB,  0.181  rbd/ws) 
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^0 

IMPACTOR 

THICKNESS 

(mm) 

X 

lARGET 

THICKNESS 

(mm) 

X/X„ 

AT 

e 

PULSE 

WIDTH* 

(ys) 

- 1 

V 

s 

SPALL 

v'ELOCITY 

(mm/us) 

0.49 

1.01 

2.06 

0.21 

0.135 

1.49 

3.04 

2.04 

0.63 

0.110 

3.04 

6.07 

2.00 

1.28 

0.095 

0.40 

3.85 

9.63 

0.17 

0.180 

Spall  data  for  copper  reported  by  Smith  showed  similar 
time  dependence  as  the  present  work  but  a  lower  incipient 
spall  level.  This  is  probably  due  to  a  softer  material 
(Rp  =  68  for  Smith's  copper)  and  an  incipient  spall  criteria 
based  on  initial  void  appearance. 


_ _ _  _ _ _ 

MANUFACTURING  OCVELORMENT  •  GENERAL  MOTORS  CORPOTATION 

MSL-7  -23,  Vol.  V 

SPALL  PROFILES 

Spall  wave  profiles  were  studied  using  a  copper-into-copper 
configuration,  as  in  the  wave  profile  tests,  but  with  a  free 
rear  surface  on  the  target  for  complete  reflection  of  the  com¬ 
pressive  front.  Measured  profiles  for  unattenuated  pulses  are 
shown  in  Figure  28,  with  photomicrographs  of  the  recovered  and 
sectioned  targets.  Also  shown  for  reference  are  the  incipient 
and  complete  spall  levels  for  this  impactor/ target  combination 
as  determined  from  the  recovery  tests.  Evidence  of  fracture 
in  the  target  is  shown  as  a  reversal  of  the  measured  release 
wave.  This  occurs  as  the  rarefaction  wave  is  reflected  from 
the  spall  plane  as  a  compressive  wave  which  then  reaches  the 
rear  surface  and  gives  it  a  positive  velocity  increment.  As 
this  "trapped"  wave  reverberates  back  and  forth  between  the 
spall  plane  and  the  rear  surface,  a  sinusoidal-type  motion  is 
observed.  This  motion  eventually  damps  out  due  to  dispersion 
and  attenuation  in  the  material  and  to  the  effects  of  edge 
rarefactions . 

The  amount  of  pullback  (decrease  in  free  surface  velocity  to  the 
point  of  first  reversal)  may  be  related  to  the  tensile  strength 
at  the  spall  plane  of  the  material  for  a  given  impact  geometry. 

(49  \ 

An  empirical  relation  for  this  is  given  by: 

Cg  =  pC(Au^g/2)  (26) 

where  o  is  tensile  strength,  p  and  C  are  local 
s 

density  and  longitudinal  wave  speed,  and  Au^^ 
is  pullback. 
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The  data  in  Figure  28  give  a  =  21  kbcu:.  Note  that  the  pull- 

s 

back  is  relatively  insensitive  to  peak  compressive  stress. 

This  calculated  spall  strength  is  not  constant  for  a  given 

material  since  it  is  proportional  to  which  has  been  found 

to  change  with  such  factors  as  pulse  width,  pulse  shape  and 

(581 

target  tenperature .  ' 

The  spall  profile  for  an  attenuated  pulse,  where  the  time 
duration  of  the  peak  tensile  stress  is  expected  to  be  less 
than  for  the  unattenuated  case,  is  shown  in  Figure  29.  The 
pullback  for  this  case  is  0.16  mm/ysec  giving  a  spall  strength 
of  'v  35  kbar,  consistent  with  time-dependency  of  spall  fracture 
shown  in  Figure  27. 


TIME  (Msec) 

Figure  29  Attenuated  Spall  Wave  Profile  in  Copper 
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SUMMARY 


The  dynamic  properties  of  OFHC  copper  have  been  measured  under 
uniaxial  stress  and  uniaxial  strain  conditions.  The  material 
tested  was  1/2-hard  bar  stock  with  64  ym  grain  size  and  78Rp 
hardness.  Material  response  to  compressive  uniaxial  stress 
deformation  approximated  elastic-perfectly  plastic  behavior. 
There  was  slight  rate  sensitivity  with  2.1  kbar  yield  at 
0.005/sec  and  3.2  kbar  yield  at  950/sec.  Annealing  reduced 
the  yield  to  0.25  kbar  but  increased  the  work-hardening  ex¬ 

ponent  by  almost  an  order  of  magnitude. 

Longitudinal  and  shear  wave  velocities  were  measured  as  func¬ 
tions  of  temperature  and  pressure,  giving  the  following  (in 
mm/ysec) : 

=  4.766  -  (0.000434  X  T) 

Cg  =  2.253  -  (0.000300  X  T) 

(T  in  ®C) 

=  4.757  +  (0.0  585  X  P) 

Jj 

Cg  =  2.247  +  (0.00243  X  P) 

A 

(P  in  kbar) 

Principal  elastic  constants  were  evaluated  at  T  -  20*'C  and 
P  s  0,  including: 
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Bulk  Modulus  *  Kq  =  1418  kbar 


C 

Shear  Modulus  =  Gq  =  450  kbar 

Poisson's  Ratio  =  0.356 

Elastic  wave  velocity  data  were  also  used  in  determining  the 
Gruneisen  parameter  y  and  the  Debye  ten^erature  6: 

Y  =  2.13  @  T  =  20“C  and  P  =  0 
=  -  O.OOOTe/^C  @  T  *  20‘>C 

IJ  =  -  0.0020/kbar  @  P  =  0 
0  =  331®K 

The  shock  wave  hugoniot  was  determined  experimentally  to  be: 
o„  =  1.5  +  344  u  +  143  u^  (o„  <  200  kbar) 

n  p  p  n 

Assumption  of  an  ideal  elastic-plastic  wave  structure  with 
an  8  kbar  hugoniot  elastic  limit  and  equilibrium  initial  and 
final  states  gave: 

»  3.1  +  1308u  +  3520w^ 

H 
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P„  =  1308  y  +  3520  y^ 

n 

U-  =  3.915  +  1.53  u 
S  p 

The  hydrostat  as  determined  from  elastic  constants  data  was: 


P 


H 


1025 


1383  In 


+  1611 


(Keane  equation  of  state) 


Measurements  of  compressive  wave  development  showed  the  wave 
front  to  be  characterized  by  a  non-steady-state  elastic  portion 
and  a  transition  to  a  spreading  plastic  wave.  Elastic  precursor 
attenuation  was  rapid,  decaying  80%  to  4  kbar  after  only  2  mm 
of  travel.  Dispersion  of  the  elastic  wave  was  so  great  that 
meaningful  values  for  the  precursor  level  could  not  be  obtained 
for  propagation  dist^'nces  greater  than  "v  13  mm.  Copper  showed 
a  relatively  low  shock  wave  attenuation  rate,  with  a  41  kbar, 

0.5  ysec  pulse  attenuating  only  12%  at  X/Xq  =  11. 

The  spall  behavior  was  determined  for  elastic  pulse  widths  of 
0.21,  0.63  and  1.28  ysec.  The  incipient  spall  threshold  expressed 
as  impact  velocity  increased  with  decreasing  pulse  width.  Copper 
showed  a  very  well-defined  spall  plane  with  void  formation  and 
growth  confined  to  a  relatively  narrow  band.  For  unattenuated 
pulses,  the  spall  plane  was  located  approximately  one  impactor 
thickness  from  the  target  rear  surface.  The  fracture  surface  for 
completely  spalled  copper  showed  ductile  failure  characterized  by 
equiaxed  dimples  with  evidence  of  serpentine  glide  at  the  void 
walls . 


72 


MANUFACTURING  DEVELOPMENT  •  GENERAL  MOTORS  CORPORATION 

MSL-70-23,  Vol.  V 

ACKNOWLEDGMENTS 


The  authors  wish  to  thank  Dr.  T.  E.  Michaels  and 
Mr.  A.  R.  McMillan  for  contributions  to  the  equation 
of  state  and  wave  propagation  studies,  Mr.  F.  L.  Schierloh 
for  assistance  in  the  strain  rate  work,  and  Dr.  B.  Taylor 
for  the  transmission  electron  microscopy.  Also,  the  assist¬ 
ance  of  the  following  Materials  and  Structures  Laboratory  per¬ 
sonnel  in  the  experimental  work  is  gratefully  recognized: 

Mr.  B.  Bielski  and  Mr.  J.  Wigton  for  specimen  preparation; 

Mr.  F.  Feistmann,  Mr.  F.  Stemczynski  and  Mr.  E.  Miller  for  EOS 
and  wave  profile  tests;  Mr.  J.  Havens  for  ultrasonics  measure¬ 
ments;  Mr.  C.  Woodcock  for  spall  tests;  Mr.  J.  Bor  Ler  and 
Mr.  H.  Simmons  for  strain  rate  tests;  Mr.  C.  Whitchurch  and 
Mr.  L.  Seitz  for  metallographic  work;  and  Mrs.  W.  Trippen  and 
Mrs.  C.  Strong  for  data  analysis. 


73 


MANUFACTURING  DEVELOPMENT  •  GENERAL  MOTORS  CORPORATION 


MSL-70-23,  Vol.  V 


REFERENCES 


1.  "A  Selected  Bibliography  on  Dynamic  Properties  of  Materials", 
compiled  by  Christman,  D.  R.,  General  Motors  Corporation, 
Manufacturing  Development,  Materials  and  Structures  Labor¬ 
atory,  DASA  2511,  June,  1970  (AD  710823). 

2.  Metals  Handbook,  8th.  Ed.j  Vol.  1;  Propei“ties  and  Selection 
of  Metals,  Lyman,  T.,  Ed.,  American  Society  for  Metals, 

Metals  Park,  Ohio,  1961. 

3.  Handbook  of  Chemistry  and  Physics,  45th.  Ed.,  Weast,  R.  C., 
Ed.,  The  Chemical  Rubber  Co.,  Cleveland,  Ohic,  1964. 

4 .  Thermophysical  Properties  of  High  Temperature  Solid  Materials, 
Vol.  1:  Elements,  Touloukian,  Y.  S.,  Ed.,  The  McMillan  Co., 
New  York,  N.Y.,  1967. 

5.  Babcock,  f.  G.,  Langan,  J.  J.,  Norvey,  D.  B.,  Michaelr,  T.  E. 
and  Schierloh,  F.  L.,  "Characterization  of  Three  Aluminum 
Alloys",  General  Motors  Corporation,  Manufacturing  Develop¬ 
ment,  Materials  and  Structures  Lciboratory,  MSL-70-12,  May, 
1970. 

6.  Maiden,  C.  J.  and  Green,  S.  J.,  "Compressive  Strain-Rate 
Tests  on  Six  Selected  Materials  at  Strain  Rates  from  10“3  to 
104  In/In/Sec",  ASME  Trans.- J.  Appl.  Mech.,  Ser.  E,  Vol.  33, 
p.  496-504,  1966. 

7.  Kumar,  A.  and  Kumble,  R.  G.,  "Viscous  Drag  on  Dislocations 
at  High  strain  Rates  in  Copper",  J.  Appl.  Phys.,  Vol,  40, 
p.  3475-3480,  1969. 

8.  Green,  S.  J.,  Maiden,  C.  J.,  Babcock,  S.  G.  and 
Schierloh,  F.  L.,  "The  High  Strain- Rate  Behavior  of  Face 
Centered  Cubic  Materials",  General  Motors  Corporation, 
Manufacturing  Development,  Materials  and  Structures  Labor¬ 
atory,  MSL-69-36,  Oct.,  1969. 

9.  Wonsiewicz,  B.  C.  and  Chin,  G.  Y.,  "Plain  Strain  Compression 
of  Copper,  Cu6WtPctAl,  and  Ag4WtPctSn  Crystals",  Met.  Trans., 
Vol.  1,  p.  2715-2722,  1970. 

10.  Morrow,  J.,  "Cyclic  Plastic  Strain  Energy/  and  Fatigue  of 
Metals",  Internal  Friction,  Dang>in^  and  Cyclic  Plasticity, 

STP  378,  American  Society  for  Testing  and  Materials,  p.  45-87, 
1965. 


74 


MANUFACTURING  DEVELOPMENT  •  GENERAL  MOTORS  CORPORATION 


MSL-70-23,  Vol.  V 


11.  Hecker,  S.,  Private  Communication#  General  Motors  Research 
Laboratories,  1970. 

12.  SalEuna#  K.  and  Alers#  G.  A.#  "Third  Order  Elastic  Constants 
of  Copper  at  Low  Temperature",  Phys.  Rev.,  Vol.  161, 

p.  643-680,  1967. 

13.  Chang,  Y.  A.  and  Himmel,  L,  "Tenperature  Dependence  of  the 
Elastic  Constants  of  Cu,  Ag,  and  Au  Above  Room  Temperature", 
J.  Appl.  Phys.,  Vol.  37,  p.  3567-3572,  1966. 

14.  Kolsky,  H.,  Stress  Waves  in  Solids,  S1098,  Dover  Publica¬ 
tions,  Inc.,  New  York,  N.  Y.  1963. 

15.  Mason,  W.  P.,  "Acoustic  Properties  of  Solids",  American 
Institute  of  Physics  Handbook,  2nd.  Ed.,  p.  3-28  -  3-97, 
1963. 

16.  Barsch,  G.  R.  and  Chang,  Z.  P,,  "Adiabatic,  Isothermal, 
and  Intermediate  Pressure  Derivatives  of  the  Elastic  Con¬ 
stants  for  Cubic  Symmetry",  Phy.  Stat.  Sol  ,  Vol,  19, 

p.  129-151,  1967. 

17.  Walsh,  J.  M.,  Rice,  M.  H.,  McQueen,  P.  G.  and  Y/irger,  F.  L. , 
"Shock-Wave  Congress tons  of  Twenty-Seven  Metals,  Equation 
of  State  of  Metcils",  Phys.  Rev.,  Vol.  108,  p.  196-216,  1957. 

18.  Slater,  J.  C. ,  Introduction  to  Chemical  Physics,  McGraw-Hill 
Book  Co.,  New  York ,  N .  Y . ,  1939 . 

19.  Dugdele,  J.  S.  and  MacDonald,  D.  K.  C,,  "The  Thermal  Expan¬ 
sion  of  Solids",  Phys.  Rev.,  Vol.  89,  p.  832-834,  1953. 

20.  Schreiber,  E.  and  Anderson,  0.  L.,  "Pressure  Derivatives  of 
the  Sound  Velocities  of  Polyccystailine  Alumina",  Am.  C^rara. 
Soc.  J.,  Vol.  49,  p.  184-190,  1966. 

21.  Anders.^n,  0.  L.  and  Dienes,  G.  J.,  Non-Crystalline  Solids 
(Chapter  18),  Ed.  by  Frechette,  V.  D- ,  John  wiiley  anci 
Sons,  Hew  York,  1960. 

22.  Schuele,  D.  E.  and  Smith,  C.  S.,  "Low  Temperature  Thermal 
Expansion  of  Rbl",  J.  Phys.  Chem.  Solids,  Vol.  25,  p.  801- 
814,  1964. 

23.  McQueen,  R.  G.  and  Harsh,  S.  P.,  "Equation  of  State  for 
Nineteen  Metallic  Elements  fxcmi  Shock-Wave  Measurements 
to  Two  Megabars",  J.  Appl.  Phys.,  Vol.  31,  p.  1253-1269 . 
1960. 


75 


manufacturing  development  •  GENERAL  MOTORS  CORPORATION 


MSL-70-23,  Vol.  V 


24.  Anderson,  0.  L.,  "A  Sin^lified  Method  for  Calculating  the 
Debye  Temperature  from  Elastic  Constemts",  J.  Phys.  Chem. 
Solids,  Vol.  24,  p.  909-917,  1963. 

25.  Alers,  G.  A.,  "Use  of  Sound  Velocity  Measurements  in  Deter¬ 
mining  the  Debye  Temperature  of  Solids",  Physical  Acoustics, 
Vol.  Ill,  Part  B  (Lattice  Dynamics),  Academic  Press,  New 
York,  N.y.,  p.  1-42,  1965. 

26.  Munson,  D.  E.  and  Barker,  L.  M.,  "Dynamically  Determined 
Pressure-Volume,  Relationships  for  Aluminum,  Copper  and 
Lead",  J.  Appl.  Phys.,  Vol.  37,  p.  1652-1660,  1966. 

27.  Isbell,  W.  M.,  Shipman,  F.  H.  and  Jones,  k.  H.,  "Hugoniot 
Equation  of  State  Measurements  for  Eleven  Materials  to 
Five  Megabars",  General  Motors  Corporation,  Manufacturing 
Development,  Materials  and  Structures  Laboratory,  MSL-68-13, 
December,  1968;  also  see  MSL-71-12,  April,  1971. 

28.  Overton,  W.  C.,  "Relation  Between  Ultrasonically  Measured 
Properties  and  the  Coefficients  in  the  Solid  Equation  of 
State",  J.  Chem.  Phys.,  Vol.  37,  p.  116-119,  1962. 

29.  Anderson,  0.  L. ,  "The  Use  of  Ultrasonic  Measurements  Under 
Modest  Pressures  to  Estimate  Compression  at  High  Pressures", 
J.  Phys.  Chem.  Solids,  Vol.  27,  p.  547-565,  1966. 

30.  Lamberson,  D.  L.,  "The  High  Pressure  Equation  of  State  of 
Tantalum,  Polystyrene  and  Carbon  Phenolic  Determined  from 
Ultrasonic  Velocities",  Doctoral  Dissertation,  Air  Force 
Institute  of  Technclogy  (AFIT-SE) ,  Wright-Patterson  AFB, 

Ohio,  Marc^ ,  1969. 

31.  MacDonald,  J.  R.,  "Review  of  Some  Experimental  and  Analytical 
Equations  of  State",  Rev.  Mod.  Phys.,  Vol.  41,  p.  316-349. 
1969. 

32.  Birch,  F.,  "The  Effect  of  Pressure  Upon  the  Elastic  Para¬ 
meters  of  Isotropic  Solids,  According  to  Hurnaghan's  Theory 
of  Finite  Strain",  J.  Appl.  Phys.,  Vol.  9,  p.  279-288,  1938. 

33.  Birci: .  F.,  "Elasticity  and  Constitution  of  the  Earth's  In¬ 
terior",  J.  Geophys.  Res.,  Vol.  57,  p.  227-286,  1952. 

34.  Murnaghan,  F.  D.,  "The  Comjressibility  of  Media  Under  Ex¬ 
treme  Pressures",  Proc.  Nat.  Ac.  Sci . ,  Vol.  30.  p.  244-247. 
1944. 


76 


manufacturing  OEVELORMINT  •  GENERAL  MOTORS  CORPORATION 


MSL-70-23,  Vol.  V 


35.  KecUie/  A.,  Australian  J .  Phys . ,  Vol.  7/  p.  323,  1954. 

36.  Anderson,  0.  L.,  "On  the  Use  of  Ultrasonic  and  Shock-Wave 
Data  to  Estimate  Compressions  at  Extremely  High  Pressures", 
Phys.  Earth  Planet.  Intorio..8  ,  Vol.  1,  p.  169-176,  1968. 

37.  Duvall,  G.  E.  and  Zwolinski,  B.  J.,  "Entrof ic  Equations  of 
State  and  Their  Applications  to  Shock  Wave  Phenomena  in 
Solids",  Acoust.  Soc.  Am.  J.,  Vol.  27,  p.  lOS-^-lCSB,  1955. 

38.  Zel'dovich,  Y.  B.  and  Raizer,  Y.  P.,  Physj.cs  of  Shock  Waves 
and  High-Temperature  Hydrodynamic  Phenomena,  Vol.  II,  Ac- 
ademic  Press,  New  York,  N.Y.,  p.  688-709,  1967. 

39.  Ruoff,  A.  L.,  "Linear  Shock-Velocily-Particle-Velocity  Re¬ 
lationship",  J.  Appl.  Phys.,  Vol.  38,  ^  4976-4980,  1967. 

40.  Pastine,  D.  J.  and  Piacesi,  D.,  "The  Existence  and  Implica¬ 
tions  of  Curvature  in  the  Relation  Between  Shock  and  Par¬ 
ticle  Velocities  for  Metals’',  J.  Phys.  Chem.  Solids,  Vol.  27, 
p.  1783-1792,  1966. 

41.  Smith,  C.  S.,  "Metallographic  Studies  of  Metals  after  Ex¬ 
plosive  Shock",  AIME  Trans ♦ ,  Vol,  212,  p.  574-589,  1958. 

42.  Appleton,  A.  S.  and  Waddington,  J.  S.,  "The  Importance  of 
Shock  Wave  Profile  in  Explosive  Loading  Experiments",  Acta 
Met. ,  Vol.  12,  p.  956-957,  1964. 

43.  Higgins,  G.  T. ,  "The  Structure  and  Annealing  Behavior  of 
Shock-Loaded,  Cube-Oriented  Copper",  Metall.  Trans.,  Vol.  2, 
p.  1277-1282,  1971. 

14.  Johari,  0.  and  Thcwias,  G.,  "Substructures  in  Explosively 

Deformed  Cu  and  Cu-Al  Alloys",  Acta  Met.,  Voi.  12,  p.  1153- 
1159,  1964. 

45.  Brillhart,  D.  C.,  DeAngelis,  R.  J.,  Preban,  A.  G., 

Co  3n,  J.  B.  and  Gordon,  P.,  "Quantitative  Study  of  the 
Substructure  and  Properties  of  Shock-Loaded  Copper", 

AiHE  Trans . ,  Vol.  239,  p.  836-843,  1967. 

46.  Nurr,  L.  E.  and  Grace,  F.  I.,  "Defect  Microstructure  and 
Mechanical  Properties  in  Shock-Hardened  Metals",  Exp.  Mech., 
Vol.  9,  p.  145-155,  1969. 

4/.  Jones,  O.  E.  and  Mote,  J.  D.,  "Shock-Induced  Dynamic  Yield¬ 
ing  in  Copper  Single  Crystals",  J.  Appl.  Phys.,  Vol.  4L, 
p.  4920-4923,  1969. 


77 


MANUFACTURING  DEVELOPMENT  •  GENERAL  MOTORS  CORPORATION 


MSL-70-23,  Vol.  V 


48.  Christman,  D.  R.  and  Froula,  N.  H.,  "Dynamic  Properties  of 
High-Purity  Beryllium",  AIAA  J . ,  Vol.  8,  p.  477-482,  1970. 

49.  Taylor,  J.  W. ,  "Stress  Wave  Profiles  in  Several  Metals", 
Dislocation  Dynamics,  McGraw-Hill  Book  Co,,  New  York,  N.Y., 
p.  573-589,  1968. 

50.  Taylor,  J.  W.,  "Dislocation  Dynamics  and  Dynamic  Yielding”, 
J.  Appi -  Phys.,  Vol.  36,  p.  3146-3150,  1965. 

51.  Gilman,  J.  J.,  "Dislocation  Dynamics  and  the  Response  of 
Materials  to  Impact",  Appl.  Mech.  Rev.,  Vol.  21,  p,  767- 
783,  1968. 

52.  Johnson,  J.  N.,  "Constitutive  Relation  for  Rate-Dependent 
Plastic  Flow  in  Polycrystalline  Metals",  J.  Appl.  Phys., 
Vol.  40,  p.  2287-2293,  1969. 

53.  Read,  H.  E.,  Triplett,  J.  R.  and  Cecil,  R.  A.,  "Dislocation 
Dynamics  and  the  Formulation  of  Constitutive  Equations  for 
Rate  Dependent  Plastic  Flow  in  Metals",  Systems,  Science, 
and  Software,  DASA  2638,  December,  1970. 

54.  Isbell,  W.  M. ,  Christman,  D.  R. ,  Babcock,  S.  G., 

McMillan,  A.  R.  and  Green,  S.  G. ,  "Measurements  of  Dynamic 
Properties  of  Materials,  Vol.  Ill:  6061-T6  Aluminum", 
General  Motors  Corporation,  Manufacturing  Development, 
Materials  and  Structures  Laboratory,  DASA  2501-3,  1971. 

55.  Christman,  D.  R. ,  Michaels,  T.  E,,  Isbell,  W.  M.  and 
Babcock,  S.  G.,  "Measurements  of  Dynamic  Properties  of 
Materials,  Vol.  IV:  Alpha  Titanium",  General  Motors 
Corporation,  Manufacturing  Development,  Materials  and 
Structures  Laboratory,  DASA  2501-4,  1971. 

56.  Isbell,  W.  M. ,  Christman,  D.  R.  and  Babcock,  S.  G.,  "Mea¬ 
surements  of  Dynamic  Properties  of  Materials,  Vol.  VI: 
Tantalum",  General  Motors  Corporation,  Manufacturing 
Development,  Materials  and  Structures  Laboratory, 

DASA  2501-6,  1971. 

57.  Smith,  J.  H.,  "Three  Low-Pressure  Spall  Thresholds  in 
Copper",  Dynamic  Behavior  of  Materials,  American  Scoiety 
for  Testing  Materials,  Philadelphia,  P.  264-282,  1963. 

58.  Isbell,  W.  M.  and  Christman,  D.  R.,  "Shock  Propagation 
and  Fracture  in  6061-T6  Aluminum  from  wave  Profile  Mea¬ 
surements",  General  Motors  Corporation,  Manufacturing 
Development,  Materials  and  Structures  Laboratory, 

DASA  2419,  April,  1970  (AD705536) . 


78 


MANUFACTURING  DEVELOPMENT  •  GENERAL  MOTORS  CORPORATION 


MSL-70-23,  Vol.  V 


APPENDIX  A 

EIASTIC  CONSTANTS  EQUATIONS 


Following  are  the  equations  used  in  calculating  various  elastic 

constants  and  the  adiabatic  and  isothermal  pressure  and  tempera- 

* 

ture  derivatives  of  the  bulk  and  shear  moduli!  (K  and  G)  : 

(Bulk  Wave  Velocity) 

Cj,  =V  2(l+v) 

(Sound  wave  velocity,  longitudinal 
wave  of  infinite  wave  length  in 
bar,  lateral  displacement  allowed) 

Cj^;  kj  -  8kJ  +  (24-l6a^^)k^^  +  (16a^^-16)  =  0 

where  kj  =  “l^  "  ^ 

(Rayleigh  or  Surface  Wave  Velocity) 


The  expressions  for  the  pressure  derivatives  are  simplified 
in  these  equations,  e.g.. 
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v=  [0.5-(Cg/Cj^)^]/[1.0-{Cg/Cj^)^] 
(Poisson's  Ratio) 

K  =  pCg  =  p(C^-4/3Cg) 

(Bulk  Modulus) 

G  =  PC| 

(Shear  or  Rigidity  Modulus) 

E  =  2p(l+v)  Cg 

(Young's  or  Elastic  Modulus) 


A  =  2pvCg/(l-2v) 

(Lame's  Parameter) 
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APPENDIX  B 

GRUl'IEISEN  PARAMETER  ESTIMATES 


1.  The  Griineisen  parameter  is  defined  thermodynamically  as 


For  an  isotropic  elastic  solid,  this  becomes 


Substituting  appropriate  values  for  copper  gives 
Yq  =  2.13 

2.  Slater  developed  a  relation  that  was  an  extension  of  the 
Debye  Theory  for  a  three-dimensional  continuum,  and 
assumed  Poisson's  ratio  was  independent  of  volume  and 
Y  is  the  same  for  all  vibrational  frequencies: 


which  simplifies  at  zero  pressure  to 
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where  S  is  from  U_  =  C_  +  Su  . 

O  O  p 

Ultrasonic  data  gives  S  =  1.60  and  y  =  2.53. 
Shock  wave  data  gives  S  =  1.53  and  y  =  2.39. 


3.  Dugdale  and  MacDonald  give  an  expression  for  plane  wave 
velocity  vs.  volume  in  a  three-dimensional  continuum  that 
leads  to  an  expression  for  y  differing  from  Slater's  by 
a  constant  (1/3) : 


which  simplifies  at  zero  pressure  to 


From  ultrasonic  data,  y  =  2.20. 
From  shock  wave  data,  y  »  2.06. 


4.  Anderson  and  Dienes  give  expressions  for  y  in  terms  of  the 
longitudinal  and  shear  modes  of  wave  propagation: 
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Ultrasonic  data  gives  ~  2.31  and  Vg  =  1.82. 

5.  Schuele  and  Smith  give  expressions  for  the  "low-ten^era- 
ture"  and  "high-temperature"  Gruneisen  ratio,  where  low- 
temperature  is  the  range  in  which  it  is  assumed  the  solid 
behaves  as  a  continuum  and  high-temperature  is  the  range 
where  all  vibrational  modes  are  excited  and  y  is  constant: 


Using  ultrasonic  data  and  and  y^  from  above  gives 
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6.  Perry  calculated  the  Grueisen  parameter  from  laser  inter¬ 
ferometer  data  obtained  for  electron  beam  energy  deposi- 
* 

tion.  For  polycrystalline  copper,  he  obtained  y  =  1.63 
assuming  instantaneous  deposition  and  y  -  1.93  assuming 
50  nsec  deposition  time. 

(All  of  the  ultrasonics  data  used  in  this  appendix  are  for 
zero  pressure  and  20 °C) . 


Perry,  F.  C. ,  "Thermoelastic  Response  of  Polycrystaliine 
Metals  to  Relativistic  Electron  Beam  Absorption" ,  J,  Appl. 
Phys . ,  Vol.  41,  p.  5017-5022,  1970, 


86 


M  A  N  U  F  A  (.  T  U  R  I  N  e  OfVELOPMENT  •  GENERAL  MOTORS  CORPORATION 

MSL-70-23,  vol.  V 

APPENDIX  C 

6RUNEISEN  PARAMETER,  PRESSURE  DERIVATIVE 

1.  Thermodynamic  Relation  - 

Y  -  Mi 

pCp 

-  L.  / 

Cp  \  3P  ;  T 

* 

This  can  be  sin^lified  by  using: 


/  38  \  ^  1 

^  \  ' 

1 

^  lip  1  T  1 

(  Tp^)  t  "  ■  I  (n)  p  * 

See,  e.g.  Thurston,  R.  N.,  "Ultrasonic  Data  and  the  Thermo¬ 
dynamics  of  Solids",  IEEE  Proc. ,  Vol.  53,  p.  1320-1336,  1965. 
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T 

and  K 


*0  ^  ‘'St 


P 


gives : 


At  P  =  0,  ti-is  reduces  to 
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APPENDIX  D 

DEVIATORIC  STRESS  ESTIMATES 


The  shock  wave  equation  of  state  for  copper  is  given  by: 

a„  =  1.5  +  344  +  143  u  ^ 

H  P  P  \ 

Ojj  =  3.1  +  1308y  +  3520y^  J  (D-1) 

Pjj  =  1308u  +  3520y^  (Hydrostat)  ^ 

where  y  =  (v^/v)  -  1. 


The  deviatoric  stress  (o^)  can  be  defined  as  the  difference 

between  the  dynamic  stress-compression  path  (a„  -  y)  and  the 

n 

hydrostat,  and  is  a  measure  of  the  material's  ability  to 

support  shear.  If  the  material  strain  hardens,  the  deviatoric 

stress  will  be  a  function  of  y.  In  Equation  D-1,  was 

assumed  constant  at  3.1  kbar  which  gives  P„  =  o„  -  3.1. 

n  n 

The  deviatoric  stress  can  also  be  estimated  from  uniaxial 
stress  data,  shock  wave  data  or  ultrasonics  data  as  follows: 


where  Y  is  flow  stress  under  uniaxial 
stress  conditions. 
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A.  The  flow  stress  at  low  strain  rate  from  the  uniaxial 
stress  tests  (see  Figure  3)  is  ''*3  kbar,  which  gives 
Op  =  2  kbar. 

B.  For  a  hugoniot  elastic  limit  (a  )  of  8  kbar,  the  rela- 

tion  Y  -  a  (l-2v/l“v)  gives  Y  =  3,6  kbar.  From  Equa- 
0 

tion  D-2,  Op  =  2.4  kbar. 

C.  Since  Op  can  also  be  taken  as  the  difference  between 

a„  and  ?„  at  a  compression  corresponding  to  the  hugoniot 
11  H 

elastic  limit,  one  has 


For  a  =  8  kbar:  u  =  0,004  and  P  =  5,7  kbar.  This  gives 
0  0  0 

Cjj  =  2.3  kbar. 

The  agreement  of  these  estimates  with  the  value  of  3.1  kbar 
from  Equation  D-1  is  considered  good,  taking  into  account 
the  assumptions  underlying  each  method  and  the  experimental 
inaccuracies . 

Also,  it  should  be  noted  that  the  stress-axis  intercept  is 
larger  in  the  -  y  plane  than  the  -  u^  plane.  Assuming 
linear  elastic  and  hydrostatic  behavior  up  to  the  elastic 
limit,  then: 
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0„  -  plane  - 
H  p 


u 


Ojj  -  u  plane  - 


Using  the  approximations  u  =  y  /C  and  y  z  y  /(1+y  ), 

G  G  G  G  G 

and  subtracting: 


'’0<=B‘'e<=L 


S' 


For  Og  =  8  kbar,  y^  =  0.004  and  AOj^  =1.1  kbar,  compared  to 

Ao„  =  1.6  kbar  from  a„  -  u  and  0„  -  y  given  in  Equation  D-1. 
D  n  p  n 
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APPENDIX  E 


EQUATIONS  OF  STATE,  ULTRASONICS 


Numerous  analytical  and  empirical  relations  have  been  developed 
for  determining  equations  of  state  from  ultrasonics  data.y^Some 
of  these  are  presented  here.  / 

/ 

/ 

/ 

1.  The  Birch  equation  of  state  for  the  isotherm  can  be 

* 

expressed  as : 

/ 

/ 

/ 

Pt  =  I  kJ  (  1+M  17^1  (  l-p)  -  1 


12-3  K 


where  C 


Py  =  2060 


/  [  1  /  \  2/3 

jl+y  J  -  f  1+u  I  1+1.19  (l+Mj  -1 


Another  parameter  sometimes  used  with  this  equation  of  state 


r  2/3  1 2 

a  (1+w)  -  1  J 

where  „  =  irKV  -  7k’"  +  ( 

“  8  ^O^OT  '^OT  I  ^OT  I  9  J 


However,  this  term  can  usually  be  neglected,  contributing 
<0.5%  at  100  kbar  for  copper. 
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2.  The  Murnaghan  equation  of  state  is  derived  by  assuming 
the  bulk  modulus  is  linear  with  pressure  which  gives; 


Isotherm  •>  P„  = 


Pt  =  245 


H 


.6 


Adiabat  - 


■  I  (‘") 


Pg  =  262.6 


H " 


3.  Instead  of  a  linear  bulk  modulus-pressure  relation,  Keane 


assumed: 


dK^  = 


( <  -  -iIa 


+  K®  I  dP 


S'  S' 

where  K  is  the  limiting  value  of  K  at  very  high  pressure. 


This  leads  to; 


kV 

*0^0S 


H 


ks  -  -fU 


(uu) 
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(36) 

Anderson'  '  suggested  that  this  equation  be  applied  by 

g  • 

empirically  determining  .  This  can  be  done  by  evaluat¬ 
ing  the  equation  at  a  known  high  pressure  shock  wave  point, 

S ' 

i.e.,  the  equation  is  solved  for  K  at  a  known  P_,  y 
point.  Evaluating  the  high  pressure  hugoniot  given  by 

C 

Equation  14  at  Pjj  =  Kq  =  1418  kbar  gives  y  =  0.443.  This 

value  of  P„  corresponds  to  P,,  =  1371  kbar,  using  the 

“  ^  (37) 

entropy  correction  term  suggested  by  Duvall  for  the 

S ' 

Murnaghan  equation.  Solving  the  above  equation  for 
at  Pg  =  1371  kbar  and  y  =  0.443  gives 


=  559.3 


[  (-) 


651.5  In 


* 

Huang 


gives  the  upper  bound  for 


as: 


S'  8  ^ 
OS  “  I 


2.73 


This  gives: 


Pg  =»  1025 


-  1383  In 


Huang,  Y.  K. ,  "Analytical  Aproach  to  the  Shock  Compressibility 
of  18  Cubic-Lattice  Metals",  J.  Chem.  Phys.,  Vol.  53,  p.  571-575, 
1970. 
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